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The Anhysteretic Magnetization of Permanent Magnet Alloys} 


By E. P. WoHLFARTH 
Department of Mathematics, Imperial College, London 


[Received March 6, 1957] 


SUMMARY 
Calculations on the basis of the single domain treatment are reported 
of the anhysteretic magnetization, the anhysteretic remanent magnetiza- 
tion and the initial magnetization for a simple model of a permanent 
magnet alloy, and the results are compared with experimental curves for 
Alcomax ITI (Gould and McCaig 1954). 


§ 1. InTRODUCTION 

THE anhysteretic magnetization J,—J,(H) of a ferromagnetic specimen 
is that magnetization which the specimen retains after the application 
of a direct field H on which is superposed an alternating field of initial 
amplitude H’ which is then slowly reduced to zero. The function /,(H) 
is single valued and independent of the value of H’ as long as this is 
sufficiently large (see § 4). If H is the internal field, i.e. corrected for 
the external demagnetizing factor which depends on the shape of the 
specimen, then it is found experimentally that the initial slope of the 
I (H) curve is exceedingly large. The theoretical treatment given below 
and elsewhere indicates that this slope is, in fact, infinite, and where 
a finite value is observed it has usually been assumed that there is present 
in the specimen a certain amount of internal demagnetization. This 
may be due to the presence of grain boundaries, non-magnetic inclusions 
or other structural imperfections. If the effect is characterized by an 
internal demagnetizing factor NV ,, then the initial slope of the anhysteretic 
magnetization curve is given by (Néel 1943) 


P(A agar EN po! We ou 2 (1) 


Measurements of anhysteretic magnetization curves have been reported 
from time to time for a variety of materials. Kahan (1934) obtained 
results for nickel and cobalt specimens over a range of temperatures, 
and used relation (1) to determine the temperature dependence of N,. 
Similar measurements for nickel and cobalt alloys were reported by 
Kranz (1956), who related the measured values of V, directly to some 
structural properties of the specimens. The structure of permanent 
magnet alloys has a great influence on their magnetic properties, and it is 
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here that such measurements are of special interest. Experiments 
on cobalt steel and an Fe—Ni-—Al alloy have been reported by Forrer 
(1943) and Néel et al. (1943), and on Magnico by Bulgakov (1950). These 
authors do not give full anhysteretic magnetization curves but only 
report the values of N,. Recently, however, complete curves for 
Alcomax III have been published by Gould and McCaig (1954). Apart 
from curves giving /,(H) parallel and perpendicular to the preferred 
direction of this anisotropic material, these authors also give corresponding 
results for the anhysteretic remanent magnetization J,,(H). These 
are obtained by measuring the magnetization which the specimen retains 
after the application of the usual anhysteretic treatment and subsequent 
removal of the direct field H. Curves are also given which relate the 
initial magnetization with direct field H along the two directions. All 
six curves are reproduced in fig. 2. 

In the present paper theoretical curves corresponding directly to those 
measured are obtained on the basis of the single domain treatment of 
Stoner and Wohlfarth (1948). 


§ 2. A StmpLe MopEL FoR A PERMANENT MaGnet ALLOY 


The work of Nesbitt and Williams (1955) and others shows that the 
permanent magnet properties of alloys of the Alcomax type are due to 
the shape anisotropy of single domain cobalt rich precipitates, and that 
the anisotropy resulting from the magnetic treatment is due to the 
preferential nucleation and growth of the precipitates parallel to the 
field direction (for a review see Kittel and Galt 1956, Wohlfarth 1956). 
Without entering into the finer details of the structure of these materials, 
it is possible to describe their magnetic properties by means of the following 
model, which is both physically reasonable and amenable to simple 
calculations. 

Consider an assembly of prolate spheroidal single domain precipitate 
particles each of volume Vo, saturation magnetization Ig) and shape 
anisotropy coefficient NV (m), where m is the dimensional ratio of the 
particles. Further, consider the particles to be oriented at a suitably 
defined average angle 6, with the preferred direction and 6, with the 
perpendicular direction. Ideally, 60° and 6,=—90°, but this state of 
affairs, which would lead to rectangular hysteresis cycles, is never quite 
realized in practice ; in general 6, is some small angle which may actually 
be estimated from the measurements under discussion. If the packing 
factor of the particles is denoted by p, where p=nV, with n the number 
of particles per unit volume of the material (precipitates plus matrix), 
then the saturation magnetization of the material is 7), where 


i ere ie ee). 
It is difficult to take precise account of the effects of particle interaction 
(Wohlfarth 1955 a), but in certain circumstances this may be done by 
defining an effective shape anisotropy coefficient N(m), where, formally, 


N(m)=f(p)N o(m). ye eer p48) 
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The magnetization J ee the material in an applied field H is given by 
Poel legis set ee gk ee ee CS) 
where 
H=N (M)L[yh= {f(p)/p}No(m)1 yh. a eee) 


Here j=cos ¢, where ¢ is the angle between the directions of J, and H, 
and / are the reduced magnetization and field parameters used by Stoner 
and Wohlfarth (1948). For a given value of h, 7 depends on the orienta- 
tion angle @ between H and the major axis of the particle, so that different 
I, H relations are obtained depending on whether 606, or 0). 

Since the function f(p) in relation (3) is in general difficult to determine, 
it is necessary to obtain the scale factor relating H and h in relation (5) 
in a different way. If H, is the coercive force measured parallel to the 
preferred orientation then, from (5), 


H.={f(p)/D3N(m)Io [hol - - - - - + (6) 
where /,(@) is the reduced critical field calculated by Stoner and Wohlfarth 


(table 4). Hence 
ees HN a AO iio! 2 elke weet 387) 


As shown in § 3, relations (4) and (7) may be used to dicsuss the 
anhysteretic magnetization of a material with given values of the four 
quantities J), H,, 0, 6,. The assumptions underlying the discussion are 
considered in § 4. 


§ 3. APPLICATION OF THE SINGLE DomMaAIn TREATMENT 
Figure 1 shows the reduced hysteresis cycle relating j and h for 6= 20°, 
as an example. 
Fig. | 
1.0 


, “any Lat ee 
Hysteresis cycle and initial magnetization curve for @=20°. 
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Let the reduced direct field applied to the particles be / and the reduced 
initial amplitude of the alternating field h’. If 


Whiley ey on) 


then under the influence of the alternating field the magnetization 
initially describes complete cycles, of the form shown in fig. 1. If the 
alternating field is now reduced then eventually its amplitude reaches 


a value h”, where 
h” <h+|h(0)|<h’, ib). See en 


which is insufficient for the magnetization to reverse back into a direction 
opposite to that of h. Further reversals of the alternating field with 
decreasing amplitude cause the magnetization to describe smaller and 
smaller arcs on the upper branch of the cycle, ending up at a point 7=7(/) 
on this branch when the alternating field is zero. 

The anhysteretic magnetization curve is thus made up of the vertical 
line from the origin to the remanence point j=cos 0, and that part of the 
upper branch of the cycle which lies in the first quadrant. If, after 
completion of the anhysteretic treatment, 4 is reduced to zero, the 
magnetization returns to the remanence point independently of h. ‘The 
anhysteretic remanence curve is thus made up of the same vertical line 
to the remanence point and the horizontal line passing through this 
point and lying in the first quadrant. As discussed in § 4 the two types 
of magnetization are reduced if h’, the reduced initial amplitude of the 
alternating field, is too small. 

In order to transform the reduced curves just discussed so as to be 
applicable directly to Alcomax III it is necessary to know the values of 
the four quantities J, H,, 6, 0,, as described in § 2. These are obtained 
as follows : 

Gould and McCaig’s experimental curves are shown in fig. 2. From 
curve | and the above discussion it follows that 47/)=14-1 x 10? gauss. 
The coercive force of Aleomax III is given by Anon. (1952) as 670 oersted. 
From curves 3 and 4 of fig. 2 the saturation values of the anhysteretic 
remanence curve parallel and perpendicular to the preferred direction 
are 47/,—13-0 and 4-7 x 103 gauss respectively. Hence, from the above 
discussion 


0=cos-+(13-0/14-1) 23°, 0,=cos1(4-7/14-1) 70°. 


In the present calculations the angles have been taken to equal 20° 
and 70°, for which full theoretical data are given in tables 3 and 4 of 
Stoner and Wohlfarth (1948). 

In order to take account of internal demagnetization (cf. $1) it is 
necessary to shear the calculated curves according to the value of Ne 
the internal demagnetizing factor. This is given by Gould and McCaig 
as 1/N,—103, as read off from curve | of fig. 2. The final calculated 
curves are shown in fig. 3. The overall agreement between the two sets 
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of curves of figs. 2 and 3 is good. There are some obvious discrepancies, 
and the reasons for these are discussed in § 4. 


Fig. 2 
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Experimental magnetization curves, after Gould and McCaig (1954, 
with acknowledgments to the authors and the Physical Society). 
Anhysteretic curve parallel to preferred. direction. 

Anhysteretic curve perpendicular to preferred direction. 
Anhysteretic remanence curve parallel to preferred direction. 
Anhysteretic remanence curve perpendicular to preferred direction. 
Initial magnetization curve parallel to preferred direction. 

Initial magnetization curve perpendicular to preferred direction. 


Fig. 3 
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Calculated magnetization curves ; see caption to fig. 2. 
§ 4, Discussion 
The assumptions underlying the present calculations include the 


following : 
(1) The precipitate particles are single domains. The electron 


microscope studies of Heidenreich and Nesbitt (1952) and others show that 
this assumption is reasonable. 
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(2) All the particles have the same shape. Shape variations may 
be quite large but, without fuller information, it is impossible to take 
adequate account of them. As discussed by Wohlfarth (1954) shape 
variations lead to a smoothing out of the magnetization curves, which 
would bring the calculated curves of fig. 3 into better agreement with 
those observed (fig. 2). 


(3) The orientation angles of the particles are fixed at 6, and Dy. 
Fluctuations about these two orientations may well occur; these again 
would cause the calculated curves to be smoothed out. 


(4) The matrix of the permanent magnet material is non-magnetic. 
It would be difficult to take account of any weak ferromagnetism of 
the matrix whose magnetization presumably changes by boundary 
movements. 


(5) The internal demagnetizing factor is isotropic. Gould and McCaig 
actually discuss the possibility of V, being anisotropic, but no conclusions 
on this point can be drawn at present since the precise origin of the 
internal demagnetization is unknown. 


(6) The initial amplitude H’ of the alternating field is large enough for 
complete cycles to be described when this field is first applied. Gould 
and McCaig observed that if H’ is too small the magnetization drops 
steeply. This decrease may be explained as follows: For an assembly 
of particles with orientation angle @ and corresponding reduced critical 
field h,(?) complete cycles are described initially only if relation (8) for 
the reduced field values is satisfied. If, however, 


h’+h<|h,()|, oe Fg ee cere Ea 
and if prior to the application of the fields the assembly is demagnetized, 
then during the process the magnetization of half the particles describes 
ares on the upper branch of the cycle and of the other half on the lower 
branch, no discontinuous reversals being possible in either direction. 
The final anhysteretic magnetization acquired in the reduced field h 
is now equal to the initial magnetization in this field, and no longer that 
corresponding to the magnetization at h on the upper branch of the cycle 
(see fig. 1). A comparison of curves | and 5 of figs. 2 and 3 (showing the 
fully developed anhysteretic magnetization curve and the initial magnet- 
ization curve parallel to the preferred orientation) shows that if (10) 
applies the decrease of the anhysteretic magnetization may be very large. 
The anhysteretic remanence, in particular, vanishes completely under 
these conditions. The situation is here completely analogous to that 
which obtains during the acquisition of the normal static remanent 
magnetization of a particle assembly after the application of a direct 
field (Wohlfarth 1955 b). In fact, a curve relating J,,(H) with (H#+H’) 
is, according to this discussion, similar in shape to that relating the static 
remanence /,(H) with direct field H. Both curves drop to zero values 
below the critical field. If variations of the particle shape, orientation 
angles etc. occur which tend to cause fluctuations of the critical field 
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throughout the assembly, the decrease of the magnetization does not take 
place suddenly. In such cases the decrease will first make itself felt 
when the combined fields (H+ H’) lie below the highest critical field of 
any of the particles, the decrease being most pronounced when (H-+H’) 
is somewhat above H,, the bulk coercivity (Wohlfarth 1955 b). Gould 
and McCaig’s reported decrease occurs when H’ is rather below 1000 oersted 
compared with H,—670 oersted. 


§ 5. CONCLUSION 


The calculations described show that the single domain treatment 
gives an adequate account of the anhysteretic magnetization of permanent 
magnet alloys. The measured curves shown in fig. 2 provide considerable 
potential information about the structure of such an alloy, Aleomax III, 
such as the degree of particle alignment (6, and 6,), fluctuations of the 
critical fields (smoothness of the curves compared with those calculated ; 
fig. 3) and internal demagnetization (N,). It may, therefore, be 
reasonable to suggest an extension of the measurements to other alloys 
with a view to correlating their structural and magnetic properties. 
In particular, information is required which might lead eventually to 
a precise determination of the origin of the internal demagnetization 
effect. 
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Magnetic Saturation in Alloys of Neodymium 
at Low Temperatures} 


By J. M. Lock 
Royal Radar Establishment, Gt. Malvern, Wores. 


[Received March 13, 1957] 


ABSTRACT 


The magnetization of alloys of neodymium with lanthanum has been 
measured from room temperature down to 1:5°K. Above about 15°K the 
susceptibility is field-independent and obeys a Curie-Weiss law modified 
by a small temperature-independent paramagnetic term, y=C/(7’—@)-+-K. 
The Curie constant C gives effective atomic moments for all alloys in 
approximate agreement with the theoretical value for the Nd?* ion, 3-62 py, 
per atom. At lower temperatures the alloys show striking saturation 
effects, and it is suggested that this can be accounted for in terms of a 
model in which the magnetic interactions between neodymium atoms are 
anti-ferromagnetic for nearest neighbours but ferromagnetic for more 
distant neighbours. 


$1. INTRODUCTION 


In a previous paper (Lock, to be published, referred to hereafter as L), it 
has been reported that a number of the rare-earth elements including 
neodymium exhibit magnetic anomalies at low temperatures which appear 
to be of an anti-ferromagnetic type. In order to study the effect of 
magnetic dilution on these anomalies susceptibility measurements have 
been carried out on alloys of neodymium with lanthanum varying in 
composition from 20° to 75% by weight of the former, at temperatures 
down to 1-3°K. 


§ 2. EXPERIMENTAL DETAILS 


A Sucksmith balance apparatus previously described (L) has been used 
for the susceptibility measurements, which have mostly been carried out 
in fields up to 4780 gauss. A few measurements, using a larger electro- 
magnet, have been made up to 10 000 gauss. 

Three specimens, containing 20-5% (alloy 1), 43-5°% (alloy 2), and 
73:°3% (alloy 3) neodymium by weight, were supplied by Johnson, 
Matthey and Co., Ltd. They were prepared by melting the pure metals 
together, and a portion of each was analysed afterwards. They were 
said to contain less than 0-003 iron impurity. X-ray powder photo- 
graphs showed that their crystal structure was predominantly hexagonal 
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close-packed. For the magnetic measurements pieces weighing 0:04 to 
0-3 g were cut from the main alloy block. 

All three specimens have been investigated from room temperature 
down to about 1-3°K. Over the upper part of the range the measurements 
were made in two fields only, 4780 gauss and 2530 gauss, and the suscepti- 
bility was deduced from the slope of the straight line joining these two 
points on a magnetization curve, plotting magnetic moment against 
applied field. At the lowest temperatures where marked saturation 
effects became apparent measurements were made in several different 
fields so that a complete magnetization curve could be plotted. 


§ 3. RESULTS 


Above about 15°K all three alloys showed almost pure paramagnetic 
behaviour (susceptibility field-independent), and fig. 1 shows a plot of 
1/y against 7’ in this range. In each case the susceptibility obeys a 


Fig. 1 
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Curie-Weiss law modified by a small term corresponding to temperature- 
independent paramagnetism, similar to the results for pure neodymium 
(L). The Curie constants per gramme for the three alloys are given in the 
table, together with the effective moments per atom of neodymium. 

At sufficiently low temperatures large saturation effects occur, and fig. 2 
shows typical magnetization curves for alloy 3 at a number of different 


Fig. 2 
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temperatures, the magnetic moment being measured in Bohr magnetons 
(up) per atom of neodymium. Two magnetization curves for pure 
neodymium are also plotted for comparison. An attempt to estimate the 
true saturation moment of alloy 2 at 1-2°K was made by using higher fields, 
but even at 10000 gauss the magnetic moment was still increasing 
appreciably, and assuming an exponential approach to saturation, was 
still 10% below its saturation value, estimated as 1-2;., per atom. 
Figure 2 shows that saturation effects decrease as the temperature is 
raised, the magnetization curves becoming practically linear for alloy 3 
between 8° and 10°K. A better estimate of this temperature is obtained 
by plotting M, against temperature, where M, is the intercept on the 
axis of M of the straight line joining two points on the magnetization 
curves corresponding to standard fields, in this case 4780 and 2530 gauss, 
and is thus a measure of the non-linearity of the curve. This has been . 
done in fig. 3, and it is seen that for all three alloys M, decreases with 
increasing temperature, falling sharply over an intermediate range, and 
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then tailing off towards zero at higher temperatures. By ignoring this 
tail and extrapolating the steeply falling portion of the curve, comparative 
values of the temperature 7’, at which non-linear effects first appear in 
the magnetization curves of the alloys may be obtained, and are given in 


the table, together with rough estimates of the saturation moment at 
e2eK. 


Fig. 3 
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The specific heat of alloy 2 has also been measured (L. M. Roberts, 
private communication). After an approximation to the lattice contri- 
bution has been subtracted, for which purpose the atomic heat of pure 
lanthanum has been used (Parkinson et al. 1951), the remaining heat 
capacity shows a spread out maximum centred on 7-:5°K, giving very 
roughly an extra entropy of 1-1 cal/deg per g atom of neodymium. 


§ 4, DISCUSSION 


The large saturation effects in these alloys were at first considered to 
be evidence of ferromagnetism, 7’, being a genuine ferromagnetic Curie 
temperature. However, this hypothesis leaves two curious features un- 
explained. Firstly, the effective atomic moments of the alloys derived 
from the calculated Curie constants are in rough agreement with the 
value 3:62, per atom calculated for the trivalent neodymium ion in a 
4T oj. state (Van Vleck 1932), and if ferromagnetic interactions between 
these ions is the cause of the low-temperature behaviour a saturation 
moment of Jgj,=36/11., per atom would be expected, more than 2} times 
the greatest value actually obtained. Secondly, it is possible to plot the 
values of 7’, against the mean separation of neodymium atoms in the 


Fig. 4 
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alloys (fig. 4), and a smooth curve through these points also passes very 
near a point representing the Néel temperature of pure neodymium, 
which is believed to become anti-ferromagnetic at 7-5°K (L). It is 
difficult to believe that such a drastic change in properties could be due 
to the small change in mean atomic separation from 3-66 A (pure Nd) to 
4:09 A (73:3% Nd). : 
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An alternative explanation of these difficulties might be that the 
saturation effects are due to a kind of unbalanced anti-ferromagnetism. 
The ordered magnetic lattice of pure neodymium could acquire a resultant 
moment if the added lanthanum atoms replaced neodymium atoms whose 
moments pointed in the same direction. If nearest neighbour inter- 
actions in neodymium are anti-ferromagnetic, but more distant neighbour 
interactions are ferromagnetic such an arrangement would seem to be 
energetically favourable, since then more of the remaining distant 
neighbour neodymium moments will have ferromagnetic order, while the 
number of antiparallel nearest neighbour moments will be unaffected. 
On the other hand, if the alloys are disordered solid solutions, as seems 
probable, then the lanthanum atoms cannot take up preferential sites 
in this way. The same effect might be observed however if it were possible 
for the moment of a neodymium atom to be transferred to a neighbouring 
lanthanum atom under the influence of the local magnetic forces in the 
metal. Such a transfer might perhaps be accomplished by electronic 
transitions between the 4f and 6s or 5d bands of a type similar to that 
postulated in the transition elements between the d and s bands to account 
for their poor electrical conductivity and for the properties of their alloys 
with the noble metals (Mott and Jones 1936). 

In the concentrated neodymium alloys this model would predict a 
maximum saturation moment of Jg,(=3-:27,) per atom of lanthanum, 
and alloy 3, containing 26-7°% lanthanum by weight, has in fact a 
saturation moment of 2-94, per atom of lanthanum, thus lending some 
support to the idea. When the concentration of lanthanum is of the 
order of 50% or greater, however, this model will no longer apply, since 
the original anti-ferromagnetic lattice of neodymium will then have 
fewer sites occupied by magnetic than by non-magnetic atoms. Under 
these circumstances ferromagnetic interactions between next-to-nearest 
and more distant neighbours should predominate, so that a saturation 
moment will still result, but with its value reduced to less than Jg,,, per 
atom of neodymium by the anti-ferromagnetic effect of the remaining 
nearest neighbours, in qualitative agreement with the observations. 

The anomalous specific heat of alloy 2 is in qualitative agreement with 
the appearance of magnetic ordering in the lattice, though the rather 
broad maximum found suggests that its onset is by no means sharp. 


§ 5, CONCLUSION 


The observed saturation effects in alloys of neodymium with lanthanum 
seem to be explicable on the assumption that the interactions between 
neodymium atoms are anti-ferromagnetic for nearest neighbours but 
ferromagnetic for more distant neighbours, provided that transfer of the 
neodymium moment to lanthanum atoms is possible to enable the 
missing moments to line up parallel. Such a change of sign of the inter- 
action between neodymium atoms is in qualitative agreement with 
theoretical work on the variation of the exchange integral with separation 
of the magnetic atoms (Bethe 1933). 
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This work should be regarded as a preliminary investigation, and it is 
clear that further experiments on alloys of varying composition are 
desirable. 
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SUMMARY 

It is the purpose of the present paper to discuss systematic methods 
of evaluating and checking certain lattice constants which arise in 
cooperative phenomena associated with crystal lattices. Many of these 
constants have been evaluated previously, but little systematic description 
of methods used exists in the literature, and there are discrepancies 
between results obtained by different authors. Several new suggestions 
are put forward, the most important being the use of a ‘no immediate 
return ’ random walk problem, and the use of finite close-packed clusters 
for checking the evaluations. 


§ 1. INTRODUCTION 


THE Mayer theory of gas condensation (Mayer and Mayer 1940, Chaps. 13, 
14) represented a ‘great step forward in the attempt of statistical mechanics 
to account for phase transitions. It was able to account for condensation 
at sufficiently low temperatures, and to indicate how a critical point arose 
tin liquid—vapour equilibrium. However, when it came to detailed evalua- 
tion of cluster integrals little progress was made beyond the first few terms 
and the problems to be faced seem indeed very formidable. Even the class- 
ification of the different types of cluster integral which arise presents great 
difficulties (Riddell and Uhlenbeck 1953), and for the sixth term, for 
example, Riddell lists 56 possible ‘ stars ’ (Riddell 1951). 

With phase transitions arising from lattice structures more progress 
is possible. In two dimensions an exact solution is available for a number 
of lattices for the Ising model in the absence of a magnetic field, and a 
number of additional features of the model, such as the spontaneous 
magnetization, can be evaluated exactly (Newell and Montroll 1953). 

The Ising model is formally equivalent to a regular binary solution 
in which the constituent atoms occupy sites on a lattice. If one of the 
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constituents is taken to be a hole we have a very crude model of 
liquid-vapour equilibrium which has been discussed by a number of 
authors, and is usually referred to as a ‘lattice gas’, From a study of 
the properties of this model it is possible to understand much more 
clearly how the phenomena associated with the critical point arise. 
By comparing this model with the Mayer theory, Yang and Lee (1952) 
were able to analyse critically the conclusions of the theory, and to point 
out where they were valid, and where they broke down. 

For the Ising model the cluster integrals are replaced by cluster sums. 
(Fuchs 1942). If interactions of all neighbours were taken into account, 
these sums would be almost as difficult to classify and evaluate as cluster 
integrals ; the simplification in the lattice model arises largely since it 
is a reasonable physical approximation to take account only of nearest 
neighbour interactions. For all regular lattices this automatically 
eliminates many types of cluster, and this is particularly true of loose 
packed lattices. The basic problem then arising is to count the number 
of ways in which particular types of configuration can occur on a given 
lattice. 

This latter problem is still one of considerable difficulty, and although 
many different authors have dealt with it, little attempt seems to have 
been made to devise systematic methods of enumerating and checking. 
As a consequence of this, many published results contain errors, and for 
several important lattices, such as the close-packed cubic, calculations 
have not been pushed very far. 

For the Ising model the cluster sums can be transformed to yield 
series expansions at low and high temperatures. The evaluation of the 
cluster sums involves * spin counts ’, i.e. a classification and determination 
of the number of configurations on a given lattice involving a finite number * 
of overturned spins. The high temperature expansions involve ‘ line 
counts ’, i.e. a classification and determination of the number of 
configurations involving a finite number of lattice bonds. 

In the present paper we shall be concerned with high temperature 
configurations ; these configurations are also of importance for high 
temperature expansions of more generalized models, such as the Ising 
model with higher values of spin, and the Heisenberg model. In particular 
the zero-field solution, and the magnetic susceptibility expansion in 
zero-field can be expressed in terms of certain irreducible constants, and 
it is our task to show how these constants can be determined with 
reliability for various lattices. We shall also demonstrate how certain 
finite symmetrical clusters can be used to verify our formulae. 


§ 2. THE High TEMPERATURE No-Fretp Expansion. 
WAKEFIELD’S Meruop 


If we consider the complete high temperature expansion for the simple 


Ising model in the presence of a magnetic field, the number of high 
temperature configurations increases very rapidly with the number of 
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lines. For 1, 2, 3...6 lines there are 1, 3, 8, 23, 66, 212 distinct 
configurations. In the absence of a magnetic field very few of these 
retain a non-zero contribution to the high temperature expansion. If a 
suitable transformation of variable is employed the partition function 
can be expressed in the following form (Newell and Montroll 1953) 


Zy=2*(cosh Ky] 14 Snir) tanh K ial} 


Here K=J/2k7', J is the interaction energy of parallel spins, q is the 
coordination number of the lattice, and n(r) is the number of closed 
graphs of r lines which can be constructed on the lattice ; only graphs in 
which an even number of lines meet at each vertex make any contribution. 
The configurations now required for the first six terms are the following : 


aN a) ) s) P< VES Ae The expansion can be 


pushed considerably further in this case, and for the simple cubic lattice 
Wakefield (1951) derived a series up to the 12th term. 


Fig. 1 


As a simple illustration of the problem of polygon counting let us first 
consider the number of quadrilaterals on the simple cubic lattice. It is 
clear that such a structure admits only of quadrilaterals that are plane 
squares. Through each site of the lattice there pass twelve squares, 
four in each coordinate plane. We neglect surface effects, assuming that, 
boundaries of the lattice are joined cyclically. The number of distinct 
squares for a lattice of NV points is thus 3N. When we come to examine 
the number of hexagons we find that there are three possible space types 
(fig. 1), and that these can be placed on the lattice in 12N, 6N, 4N ways 
respectively. This is the method of division into space types introduced. 
by Wakefield, but the process becomes very laborious after the first few 
terms ; there are 73 types of decagon, and 756 types of dodecagon. 

The chief difficulty, apart from the labour involved, is to ensure that 
no space types have been overlooked at each stage. For a lattice like 
the simple cubic it is possible, although laborious, to do this geometrically. 
We further illustrate the method by listing in Appendix I the 11 types 
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of octagon for this lattice. However for a lattice like the face-centréd 
cubic lattice the method is very cumbersome, and the polygons are more 
difficult to visualize and draw. Te 
_ Although we shall in the next section describe an alternative method 
of counting polygons, the classification of the first few polygons on a 
lattice into space types is of great use; the space-type technique is also 
the only one we possess for estimating the value of certain other 
configurations which arise in series expansions for non-zero field. 


bh 


§ 3. THe Meruop or CLosep CrRcUIT COUNTING ie 

‘An alternative method of determining the number of n-sided closed 
polygons on a lattice depends on the use of an n-step random walk on 
the lattice. We again illustrate this by the simple example of quadri- 
laterals on the simple cubic lattice. If we leave the origin along any axis, 
say Ox, which we may do in six ways, we have after this first step four 
choices for the next step that will each enable us to return to the origin 
in one non-self-intersecting path of four steps in all. Thus there are 
24N possibilities for the whole lattice, and we shall have described every 
square from each vertex in both directions or eight times. The number 
of distinct paths is thus 3. 

In general if w,, is the number of non-self-intersecting random) walks 
of n steps starting and ending at a given lattice point, then the number . 
of distinct polygons is equal to Nw,/2n (=Np,). Our problem thus 
reduces to determining u,,, but this is again a problem of classic difficulty. 
However, it is quite an easy matter to determine for various lattices the 
total number, r,,, of random walks of n steps which return to the starting 
point (see Appendix II). We have also found it possible from this to 
determine q,, the total number of random walks which return to the 
starting point, no immediate reversals being allowed (Appendix IIT) 


(this means that two successive steps —* cannot occur, although 


=<] can occur). A similar problem has recently been discussed 


by Temperley (1956) for the simple quadratic lattice, but our result is 
more general as it applies to all lattices. 
It is now easy to see that for values of n<6 only simple polygonal 


paths are possible, e.g. n=3 jae n=4 hea n=5 ey Thus 


we find that 


Iz= 6p, ] 
Gia 8p, rhs Seas, ee (2) 
q5—10p; he 


For values of n>5 self-intersecting paths will occur. Hence for n—6 
in addition to a simple hexagonal path, three paths are possible which 
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involve the description of two triangles (fig. 2). The first two correspond 
to lattice constants of 6 and 5 lines respectively, which we have denoted 
by pe. and p;,t. The simple hexagon will contribute 12p, paths to q¢5 
the last configuration, two triangles described successively, will contribute: 
6p; (the triangles described in the same sense since immediate reversals 
are not allowed). For p,, there is a choice of paths in accordance with 
the scheme shown in fig. 3.. From the vertices ABDE there will be- 


Pt AA 


Fig. 3 
A D A D 
G 
& 
B E B E 
Fig. 4 (a) 
A B 
S D 


four ways of describing the configuration (as shown or in the reverse 
direction), and from the vertex C there will be eight ways; hence the 
total contribution is 24p,,. Finally for p,;, starting from the vertex A 
or D the path can be described in two ways (fig. 4 (a)) (and reverse 
direction). Starting from B or C the path can be described in four ways 
(fig. 4 (b)) (and reverse direction). Hence we may deduce that 
Jg=12p54+ 240,.+ 12p,, + Sp3.- ie a ea 
We are thus furnished with a relation between p, and other lattice 
constants such as Peo, Psa: It is convenient to refer to these as reducible 


+ We have found it convenient to regard all closed configurations of n lines 
as lattice constants. As each new type arises we denote it by a symbol ; 
p, for a polygon, py, for any other constant of lines. All lattice constants 
are defined per site. so that the total number of distinct configurations is Vpy,. 
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lattice constants since they are obtained by putting together two lower 
order configurations with a common side or vertex ; we shall deal with 
the calculation of these reducible lattice constants subsequently. 


Fig. 4 (b) 
A B A B 
c D G D 
It will be seen that the coefficients of the various constants in the 
relation (3) can only be obtained directly after a somewhat intricate 


analysis. We shall show in the next section how finite clusters may be 
used to provide a substantial check on these coefficients. 


§ 4. Tur Use oF FINITE CLOSE-PACKED CLUSTERS 
A close-packed cluster of M spins is a set of WV spins in which every 
spin is connected with every other spin. For 4 spins this corresponds 
to a tetrahedron, for 6 spins to a hexagon with all bonds connected as 
shown in fig. 5. Close-packed clusters possess the simplifying feature 


Fig. 5 


ey 


that every spin is connected to every other in the cluster and hence 
polygonal counting problems are considerably simplified. Many of the 
problems of interest for standard lattices can be more easily considered 
for close-packed clusters, and hence a check can often be ‘obtained on 
the counting of constants for these lattices. 

A close-packed cluster of M points has coordination number (M—1) 
which we denote by (o+1) in accord with our practice for standard 
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lattices ; we shall often refer to this cluster briefly as CP(M). The 
quantities 7,, and q, defined in §3 can readily be evaluated for finite 
close-packed clusters, and the derivation is given in Appendix IV. 
Formulae (24) are polynomials in co representing the number of paths 
teturning to the origin after n steps for a close-packed cluster of (a+ 2) 
points, no immediate reversals being allowed. It is convenient to write 
these polynomials in the special form : 
In=t,(0+1)o+B,(¢+1)o(o—1)+y,(0+1o(o—1)(o—2)+.... (4) 
The number of ways per site of choosing 7 points on CP(M) is 
(o+2)(o+1)...(o—r+3)/(o4+2)r}, 

and therefore in (4) the coefficient of x, is proportional to the number 
of ways per site of choosing three points on the cluster, the coefficient 
of 8, is proportional to the number of ways per site of choosing four 
points on the cluster, the coefficient of y, is proportional to the number 
of ways of choosing five points on the cluster, and so on. 

We now refer back to formula (3) of §3 and note that it is a general 


formula applying to all lattices and to close-packed clusters. In general 
we may write 


Te AD a Sod ee aca et) 
all Deg 
ss 


The ,, (such as 12, 24, 12, 6 of (3)) may be referred to as path factors. 
The same lattice constant can occur for different values of n. Thus 
Ps, contributes 12p,, to gg, 28p,, to q,, 0 to gg and 36p,, to q. Hence 


A sq(6)—=12 Ag(6)=12 
Agi j=25 Ago (6)= 24 
(6) 
Aeloj = A3(6) =" 6 
A,,(9)=36 i 
Now in the sum represented by (5) there will occur lattice constants 
with 3, 4,5...s vertices. If we distinguish them by writing the number 
of vertices in p,, as a superscript, we have 
§ 
qd, > > Nan Dax o a S & . a (7) 
8<n V=3 


For a close-packed cluster any lattice constant with v vertices is 
proportional to the number of ways (per site) of choosing » spins in the 
cluster, i.e. 

Don Cnn (+ l)o(o—1)... (o—v+3). ee a 
Here the c¢,,” are pure numbers independent of o which can easily be 
determined by the symmetry of the configuration. Thus C4°=75 since 
a polygon can be described in either direction from any of its vertices. 
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It is not difficult to show for the remaining configurations arising in (3) 
that ¢,°=+4, Cs,4=4 and c,?=§. 
Thus for a close-packed cluster (7) becomes 


Gn= >> S doe Cine (9+ l)o(o—1) « «. (o—v-+3). ma Pe 


sn v=3 


and comparing with (4), we find that . 


3 — 
pa nee Cop Sty 
3 


Psy 
s<n 


>) Xa Cse'=Bn 
Pag" ye. b> ly. Boia EO 


s<n 


S) A ae ae vA n | 5 


giving (n—2) relations that must be satisfied by the path factors of q,,- 
For n=6 these relations suffice to determine the path factors uniquely- 
For n>6 there are not sufficient relations to determine the path factors 
in this way but we obtain (n—2) checks on the validity of the value 
determined by the direct method. We have been able to derive the 


path factors in (5) for values of n up to 9. be 


§ 5. REDUCIBLE LaTricE CONSTANTS eee. 


We use the term reducible lattice constants to denote constants which 

may be constructed by placing together two lower order configurations 
at one vertex or along one side. Examples of such reducible lattice 
constants are p;, and p,, referred to in §3; the former arises by placing 
together two triangles with a side in common and the latter by placing 
together two triangles with a vertex in common. It should be noted 
that although we have used the term reducible for these configurations, 
they may nevertheless be multiply-connected (e.g. p;,). 
- The method of calculation of such constants is best illustrated by 
examples, and we shall calculate p,,, p;, in terms of lower order constants. 
We first note that since there are Np, simple »-sided polygons on the 
lattice, and each polygon has n° vertices, np, such polygons pass through 
any site of the lattice. Similarly, if q is the coordination number of the 
lattice (=o-+1) there are Nq/2 bonds in the lattice, and each polygon 
has 7 sides; 2np,/q¢ polygons therefore pass through each bond of the 
lattice. 


We now form p,, I by putting together two triangles at a 
vertex. We can choose the first triangle in p, ways (per site), and the 


vertex in three ways. We can then choose for the second triangle any 
of the 3p, triangles passing through this vertex. In doing so we shall 
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describe p¢. twice, because of the symmetry with respect to the constituent 
triangles. But we shall also describe other configurations : thus if the 
two triangles have two vertices (that is one side) in common we shalk 
obtain pz, ; every possible p;, will be described four times, twice from 
each of the vertices at the end of the diagonal. Finally the triangles 
may have three vertices in common; we shall then describe every ps 
three times, once from each vertex.’ Hence we may write 


9p3?—= 2 ee +4D5q+3D3- hp tm le 4 CLE} 


We form p,, by placing two triangles together along a line. The first 
triangle can be chosen in p, ways (per site) and has three sides ; having 
fixed the side, the second triangle can be any of the 6p,/q¢ triangles passing 
through this side. We form p;, twice because of the symmetry of the 
configuration with respect to the triangles, and p, three times, once 
from each side. Hence we have, 


18p57/¢= 25, + 3p. Lipps MER UF 


Clearly formulae (11) and (12) suffice to determine p,, and p,;, in terms 
of p, and q. ‘ 

We now observe that relations such as (11) and (12) must also hold 
for finite close-packed clusters CP(/), for which the values of the constants 
can be much more readily evaluated (§ 4, eqn. (8)). Hence these finite 
close-packed clusters can be used as a check on equations such as (11) 
and (12). It is possible to substitute for the lattice constants the algebraic 
polynomials in o given by (8) and to check both sides algebraically. 
However, it is usually more convenient to employ numerical checks 
since if M is sufficiently large all the configurations occur on CP(M/). 
Thus for CP(5) p5,=6, Pg-=3, p3—=2 and it is easily verified that the 
relations (11) and (12) are satisfied. 


§ 6. TRREDUCIBLE LATTICE CONSTANTS 
Any constants which cannot be expressed in terms of lower order 
configurations we refer to as irreducible. The most important class of 
irreducible constants are the simple polygons. After this we must deal 
with constants involving a double bridge, such as pga, Pre> Pg. Shown in 
fig. 6. These can be counted by the space type method, but it is somewhat 
simpler to use a modified circuit counting technique. The simple polygons 


Fig. 6 


6a 
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which form the framework of the configuration must be classified into 
different space types, and for each type we must count the number of 
sites adjacent to the pair which are at the end of the double bridge. 
We illustrate by obtaining the value of yp, for the face-centred cubic 
lattice. 


Fig. 7 
A B A 
t 
D ! B 
) G Cc 
3N 12N 
(a) (b) 


I2N 6N 
(c) (d) 
(a) Flat square. (b) Flat diamond. (c) Obtuse diamond. (d) Acute diamond. 


ii 


epee this lattice there are four different types of quadrilateral which 
occur 3N, 12N, 12N, 6N times respectively (fig. 7). For type (a) there 
are four double bridges two each from AC, BD. For type (b) there are 
two double bridges from AC. For types (c) and (d) there are four double 
bridges two each from AC, BD. (This can only be seen properly with 
the aid of a model.) Hence we have a total of 108N types. But in 
carrying out this procedure we have counted each pa exactly three times 
since each of the three double links can equally well be regarded as the 
double bridge, and the remaining pair of double links as the basic quadri- 
lateral. Thus pg, is equal to 36. 
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, When all the acibls lattice constants of a given order have been 
aletenmined the analysis of ¢,, provided by an equation such as (3) enables 
us to check the determination. Thus we have been able to obtain the 
following three relations to follow (2) : 


Pe i216 —443 — 3952( 1 —3/q)] | 
ue Pa= 1al92 +793(93/9—1)(43/9—3)—74594(1—3/q)] | 
Pe+3P6a—SPge—33Pga= tel7g—594— 8095 - (13) 


—4(9,?+ 29395)(1—3/¢)+ 8457(3¢4—845)/9? 
+ 32q3(4¢3—44)/¢—493"(q—4)] (q=o+1) 


§ 7. GENERAL CONCLUSIONS 


Our main purpose in the present paper has been to indicate methods 
-of evaluating and checking lattice constant determinations which arise 
in crystal statistics. The same constants enter in a variety of problems 
connected with cooperative phenomena, and it has seemed worth while 
deyoting attention specifically to their determination. An analysis of 
‘previous results of even the most careful workers has revealed dis- 
crepancies and errors, and we have therefore concentrated particularly 
on checking our results. It is perhaps too much to expect that all 
sources of error have been removed in so complex a problem, but we 
hope that at least a substantial number have been eliminated. 

‘We hope elsewhere to publish fuller details of the evaluation of higher 
order lattice constants, and tables of the constants which we have 
-sueceeded in determining. We wish here to summarize the results of 
our investigation, and to quote some of the series expansions which have 
been derived from them. Thus for the expansion of the partition function 
-of the Ising model quoted in (1) we may write : 


n(3)= Nps ] 
n(4)=Npy | 
‘ n(3)=Npg (14) 
n(6)=3.N"ps?—N(P5q+ 2P3— Pe) 
n(7)—=N*p3p4—N (2prg+ Pop t+2P5a—P7) J 


and we have derived (8), n(9) for a general lattice. For a loose-packed 
lattice many configurations do not enter and we have obtained »(10). 
For the initial susceptibility of the Ising model at high temperatures 


we may write 


o= IE ont 1-(o— Iver" + Se" | Se dec 3) ha) 
3 


Lattice constants 


Constant aE cine saitied ae enti CP(5) CP(6) 

cubic cubic sar 

gaige 3 5 7 5 denies that + 

ee ONS 0 0 0 2 8 2 10/8 
Pat a 1 3 12 3 33 So) = 16fe 

Dst CY 0 0 0 6 168 12/5 . | 
Da J] 0 0 0 3 36 6 15 
Det = 2 22 148 15 970 0 10: 
Peart 0 0 12 0 36 2 10 
Pen oO 0 0 0 12 384 12 60: 
be P< 0 0 0 9 204 3+... 16 
Peat 0 0 0 0 2 1 a 
prt Gaia 0 0 0 42 6168 0 0 
ie — 2 Is 192 2 966 0 30, 
Be @) 0 0 0 30 2400 0 60: 
Pr ay 0 0 0 6 192 12 60 
Dra Db] 0 0 0 42 2196 0 15- 
Bie CN 0 0 0 30 2040 0 30 
Patt i: 0 0 0 0 600 0 30 
ee <> 0 0 0 0 24 2 10 
Deni A 0) O O 0 48 6 30" 


+ Irreducible. 
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where m is the magnetic moment of a single spin, y=tanh K, and the Cy 
are independent of c. We then find that 


C3=— 6p; 
Ca=— 8p, 
C5=6p3— 10p5+ 8p5, 
| Co—6p3+ 8py— 126+ 2854+ 8( Pea+ Pent Pec) 
C,= 10 ps—14p,+36p5,+16 peat 30 pep + l6p6, 
+ 8(PratPatPret PratPretPrt)+44p7¢ | 


We have also derived cg for a general lattice and c, for a loose-packed 
lattice. Expansions of this type for the Heisenberg model have been 
given by Rushbrooke and Wood (1955); the same constants enter, as 
indeed they must for all lattice models of this kind. The detailed physical 
conclusions to be drawn from these susceptibility expansions have been 
discussed elsewhere (Domb and Sykes 1957). The lattice constants 
corresponding to seven lines and less are presented in the table for a 
number of lattices. 

Finally, we quote series expansions for the partition function per spin 
for the simple cubic, body-centred cubic and face-centred cubic lattices, 
which are derived by taking the Nth root of Zy in (1) and correspond 
eae to putting N=1 in (14): 


| Simple cubic 
Z=2Acosh K)3{1+3v4+22v8+ 192v8+ 2046y+ 24 861v22, . .] 
Body-centred cubic 
Z=2(cosh K)4{14 12v4+ 148y°+ 2568v8+- 53 944y0+. . | 


Face-centred cubic 


Z=2cosh K)*[1-+ 8v3+ 33v4+ 168v5+ 96208 + 5928v7+- 38 9078 
4.268 0569+. . .] | 

The coefficient of »!° for the simple cubic lattice differs from those of 
Wakefield (1951) and Trefftz (1950) who both give 2070; the coefficient 
of v!2 differs from that of Wakefield who gives 24 943 (Rushbrooke 1952). 
However we have discussed this point privately with Wakefield who 
agrees that some of his space types were counted wrongly. The term 
in v! for the body-centred cubic lattice is new, as are the terms in v* 
and v® for the face-centred cubic lattice. The discussion of the physical 
conclusions to be drawn from the expansions (17) is deferred to another 


paper. 
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AVP PE NED lex 


Octagons on the simple cubic lattice (after Wakefield). (11 types- 
Total contribution 207N. p,=207.) 


Ck) 
6N 


AO Per HPNSE) Taxes LD 
The Number of Returns to the Origin after n Steps of a Random Walk 
This problem has been considered in detail elsewhere (e.g. Domb 1954) 
To each particular lattice corresponds a generating function (a, y, 2), 
and the number of returns to the origin after n steps is given by the 
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coefficient independent of x, y, z in [d(x, y, z)|". The generating functions 
for the simple well-known lattices and corresponding values of r, are ag 
follows : 


Simple quadratic } 
o=4(e*+a+y"1+y), 
r,=(n !)?/[(4n) !]* (nm even). 

Triangular 
$= 4x as Sas te ae | 
n= 35 y oe-g_ 19)! I 

72 0 == [(t-+q)/2 Pq! (s—q)!’ 


oat 1 51.0 (28--t)—n.; (t+q) even. 
‘Simple cubic 
o=e(z 2 +2+y '+y+21+2), 
! 
gee 4 aac ,7,8,t=0,1...»,r+s+t=v, n= 2p (even). 
Body-centred cubic 
p= a(x *+2)(y*+y)(2-*+2), 


r,=(n!)8/{(4n)!'8 (nm even). 


(18) 


Face-centred cubic 

esol) 4 ye 2) +(e (2) a) at -o)(y! +9) I, 
See et 2 ert)! (r-rt)! ee 
= 2 Tet Al idetsy PiGeto) PiGe+ey ye” et” 


Mee Looe 98: r, 8, t all odd or all even. 


T 


ee ee aor ore >—_—_— 


Ae Pa NPD xX. LT 
Random Walk with no Immediate Reversals 
We have been able to derive the following general formula for q,, in 
terms of 7,, for any lattice of coordination number (c+!) (7)=1), 


U : ie ( Del Vn Onr , 9a" a a Vn—4 
ee ee ae. 09) 


There are two independent methods of establishing this formula. The 
first and most obvious is to treat the random walk with no immediate 
reversals as a Markoff process. The generating matrix for a simple 
quadratic lattice, for example, is : 
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and gq, is the coefficient independent of x, y in spur (A”). It is easy to 
establish that formula (19) holds in this case, and the method is similar 
to that used by Temperley (1956). However, for a lattice such as ‘the 
face-centred cubic a matrix of order 12 needs to be considered, and 
although the method can be followed through it is quite involved. In 
any case since only the coordination number enters into formula (19) 
it should be possible to establish it by a method which is simultaneously 
valid for all lattices. 

An alternative approach was therefore put forward as follows. In any 
yandom walk on a lattice starting at the origin let 7,,(a, b,c) represent 
the number of paths in which the walker is at (a, b,c) after n steps, 
and let q,(a, b,c) represent the corresponding number in a walk with 
no immediate reversals. Then 7,,(a, b,c) can be divided into mutually 
exclusive classes in which the first reversal —$— occurs after 7 \steps ; 
each particular class thus consists of the Ee r steps with no reversal, 
then a reversal which can occur in o ways and then the remaining (n—r—2) 
steps of an ordinary random walk. Hence we may write : 


r (a, b, c)=q,(a, 6, c)+oq,_.(a, 6, c) +o Da Gl h ab. Cr (a aD a 


ab vor 


ae sy Gos AGH ie c')rs (a”, ine c”) 


wb’ 


AS nee > Gal @ . O46! Wee. Gearo. ee) 


abe 


+(o+l)r, (a,b,c), . . (21) 
a’ +a"=a } | | 
b’+b"=b 
OSs | 


Introducing generating functions : 
lie (ane) =>) (Gn Deyn oee 


Ot, G, 


V(t, Y,2)= SY nla, 5, c)aty bz? 
a, b,c, 
we find that (21) leads to 
B=Q,. +7{(Qn-oho +Qn shi t+. : FQ E,,3}+(oFI)R, 2 (22) 


and hence we can solve for @,, in terms of R,. Further details of the 
method (which has been developed in conjunction with M. EK. Fisher) 
will be published elsewhere. 


Ale Neb Ix ery 
Random Walk on a Close-Packed Cluster 


Select one point as origin O and suppose a random walk is started at QO. 
Then at each stage there are (+1) possibilities, and the total number 
of paths after n steps is (o+1)". Of these the number at O is by 
definition 7,,. Since the cluster is perfectly symmetric in regard to the 
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constituent points, the number of arrivals at each of the (o+1) other 
points must be 


[(o+1)"—r,]/(o+1). wes cer 8 (25) 
If now one more step is taken the number of arrivals at O will be (o+1) 
times the expression (23), and this is by definition 7, ae ~Lhus 
Tnry=(o+1)"—,7, 
and hence the 7,, can readily be determined as polynomials in o. 
The second method of Appendix III can be adapted to apply to 


close-packed clusters. Hence using formula (19) the q, can readily be 
determined as polynomials in c. The first few polynomials are as follows : 


Wm=9 

q2—=0 

Ig=o(o+1) 

qa=o(o+1)(o—1) (24) 
75—=0(o+1)(o?— 30-4 2) 

dg=a(o+ 1)(o?— 307+ 50—2) 

qg=0(o+-1)(o4— 303+ 70?— 85+ 3) 
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ABSTRACT 


Very thin layers of several metals have been deposited onto smooth 
(111) surfaces of silver single crystals. The initial stages of growth of 
these deposits were followed by electron diffraction. Lead, thallium 
and silver grew in the form of oriented monolayers, while three-dimensional 
crystallites of gold and nickel were found in the earliest detectable stages. 
Thin deposits up to about 8 A of both tin and antimony possessed an 
amorphous structure only. Spontaneous crystallization of these films 
was observed and thicker films were wholly crystalline. In several cases 
a proportion of randomly oriented deposit was formed which predominated 
when the thickness of the film was increased. This shows the importance 
of examining the initial stages of the process. No evidence for 
pseudomorphic growth was observed. 


§ 1. INTRODUCTION 


Many results are now available on the structure of thin films deposited 
on single-crystal substrates. Such deposits are often oriented relative 
to the substrate; a comprehensive review of these epitaxial growths 
has recently been given by Pashley (1956). In very few of the 
experiments, however, have deposits of less than 5 A in mean thickness 
been examined, and information concerning the earliest stages in the 
formation of an oriented overgrowth has not been obtained. In most 
cases it is not known whether the first atoms to arrive on the substrate 
spread over the surface to form an oriented monolayer, or whether they 
coalesce to form three-dimensional crystallites. It is important to obtain 
such information, so that in any theory of oriented overgrowths a model 
is used which is in agreement with the experimental observations. 

To make these observations, it is necessary to perform experiments 
under highly controlled conditions, and to record the growth of the 
deposit continuously in the range of thicknesses from 0 to 54. The 
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evaporation technique for producing the deposit on the substrate in situ 
in an electron diffraction camera is the best way to obtain a controlled 
rate of deposition, and the growth of the film can then be followed from 
changes in the diffraction pattern ; apparatus built for this purpose has 
already been described by Kehoe et al. (1954, 1955). -The growth of 
copper deposits on silver single crystals has been studied in this way, 
and radioactive copper was used to determine the thickness of the 
deposits (Newman and Pashley 1955). It was found that a surface 
film with a mean thickness equal to that of a monomolecular layer, 
deposited on a flat substrate surface, was easily detectable by electron 
diffraction. In fact the thinnest layer of copper which gave rise to a 
distinct diffraction pattern was only 0-8 Ain mean thickness. The copper 
formed small oriented crystallites about 12 4 in height and 60 A in width, 
even in the earliest stages of growth; there was no evidence for the 
formation of a monolayer as would be required in the theory of Frank 
and Van der Merwe (1949 a, b). 

The purpose of the present work is to extend these observations to the 
initial stages of growth of other metal deposits. Silver single crystals, 
prepared by evaporation onto heated mica, have again been used as the 
substrates. These specimens possess surfaces which are essentially 
smooth on an atomic scale, which gives the optimum condition for the 
observation of the deposits. In addition, no difficulty is experienced 
due to charging of the specimens in the electron beam, as in the case of 
non-conducting crystals. The growth of thin films of several metals 
deposited at a rate of 1 to 10 A4/min onto these substrates has been 
examined in detail. 


§ 2. EXPERIMENTAL TECHNIQUES 


All evaporations were carried out in the diffraction camera at a pressure 
of better than 10°*mm Hg. Most metals were evaporated from either 
tungsten or molybdenum wire spirals, although lead and thallium were 
evaporated from a nickel wire (see Caldwell 1941). At the low rates of 
deposition which were used, there was no appreciable heating of the 
silver substrates which were at room temperature in all cases. 

Several methods were used to estimate the thickness of the deposits. 
The average rate of evaporation was calculated from the time required 
to evaporate a known mass of the metal. The rate of deposition could 
only be obtained from this value if it was assumed that there was no 
re-evaporation from the surface of the silver substrates, as discussed by 
Frenkel (1923). In addition, the diffraction patterns from the deposits 
were compared with those obtained from copper deposits and silver 
bromide overgrowths in experiments where radioactive tracers were 
used to measure the average thickness directly (see Newman and Pashley 
1955). This comparison technique has one serious disadvantage. If 
the initial layers of a deposit formed a pseudomorphic overgrowth, it 
would be almost impossible to detect whether the layer had in fact formed, 
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since there would be little or no change from the original substrate 
diffraction pattern. In order to overcome this difficulty, another 
technique was employed. A cleaved rocksalt specimen was mounted 
next to the silver substrate on the specimen holder and thus deposits 
of approximately equal thickness were formed on both surfaces. The 
deposit formed on the rocksalt could be examined at any time by a simple 
deflection of the electron beam from the silver crystal. The rate of 
deposition could be calibrated in the following way. It has been found 
that a layer of copper 2 A in thickness, deposited onto rocksalt at room 
temperature, gave rise to an observable ring pattern on the fluorescent 
screen (see Newman 1955). The limiting thickness for any other metal 
to give a similar pattern was calculated from the difference in scattering 
power between the metal and copper. The rate of deposition onto the 
silver specimen in the initial stages was thus found from the time taken 
before a ring pattern became visible from the same thickness of deposit 
on the rocksalt sample. This technique was particularly valuable in the 
case of silver and gold deposits because the diffraction patterns from these 
films were coincident with that of the substrate. 

The methods of thickness estimation described above are obviously 
far less accurate than the radioactive tracer technique. The accuracy 
of the present measurements is not, however, considered to be in error 
by more than +50% 

Estimates of the size of crystallites of the deposits were made from the 
width and height of the diffraction spots and the usual reservations apply 
to the values given (see Pinsker 1953, p. 201). 


§ 3. RESULTS 
3.1. Lead 


The presence of a lead deposit was first detected in the diffraction 
pattern by the gradual appearance of continuous streaks of intensity 
inside the streaked spots due to the silver substrate. (A diffraction pattern 
from a typical silver specimen has already been published—see Newman 
and Pashley 1955.) The orientation was (111) lead parallel to (111) 
silver with [110] lead parallel to [110] silver. These effects were observed 
when the mean thickness of the deposit was less than 1 A. As the thick- 
ness was increased up to about 3 to 4 A, these streaks became much more 
intense, but retained their continuous nature perpendicular to the 
surface, even when the specimen was viewed at different angles of 
incidence in the electron beam. In addition, extra weaker streaks 
appeared due to double scattering via the silver substrate. A typical 
pattern is shown in fig. 1 (Plate), in which the streaks about the centre 
line of spots are due to the double scattering process. 

Since no intensity changes could be seen along the length of the streaks, 
it was concluded that the lead had formed a monomolecular layer on the 
silver surface. If the lateral dimensions of. this layer had been large, 
the streaks would have been exceedingly sharp, but measurements of the 
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width of the streaks corresponded to a crystal size of 80 A. The deposit 
is thus thought to be in the form of islands of monolayer. The lattice 
spacing of these deposits was determined relative to the silver base 
and a value of the closest distance of approach of a,—3-42-+0-02 A 
was obtained from measurements on seven plates. This value would 
correspond to a cube length of aj=4:84-10-03 4 which is about 2% 
smaller than the x-ray value of 4:92 4. A contraction of this order 
may be explained by the fact that, in a monolayer structure, the coordina- 
tion number is only six instead of the usual value of twelve in a 
three-dimensional lattice. There was no evidence to suggest that the 
lead had the same lattice constant as the silver in the earliest stages of 
growth of the layer. 

As the thickness of the deposit was increased beyond 3 to 4 A, variations 
in intensity appeared along the formerly continuous streaks. Gradually, 
definite spots were resolved, although they were still linked by streaks. 
The spot pattern indicated that equal amounts of lead were oriented 
with [110] and [110] directions parallel to the silver [110] direction 
(this parallel and anti-parallel orientation will hereafter be termed double 
positioning). In addition, rings due to randomly oriented lead were 
observed. When the thickness was increased to about 6 A, the pattern 
from the silver was almost completely obscured as shown in fig. 2 (Plate). 
An unusual feature of this pattern is that the regions of greatest intensity 
on the side diffraction streaks and spots occur at increasing distances 
from the shadow edge with increasing distance from centre line. A 
tentative explanation may be that the lead is not oriented strictly parallel 
to the silver but that there are two orientations such that the lead [110] 
and [110] directions make small positive and negative angles with the 
silver [110] direction. This would lead to two diffraction patterns each 
slightly off azimuth with respect to the original silver pattern which was 
set accurately before deposition was started. The composite pattern 
would then be expected to resemble that shown in fig. 2 (Plate). 

Lead was deposited further until the mean thickness of the layer was 
about 60 to 100 4. These layers gave rise to diffraction patterns which 
showed sharp spots due to lead in double positioning and a sharp set of 
rings due to randomly oriented lead. There were also extra faint spots 
due to twinning on {111} planes. The height and width of the diffraction 
spots indicated a crystallite height of about 40 A and a width of 80 A, 
while the width of the rings corresponded to crystallites of about 100 A 
in diameter. These values probably represent only lower limits due to 
the lack of penetration of the electron beam into the crystallites because 
of the high scattering power of lead. 

These results were obtained with a rate of deposition of about 10 A/min. 
If a rate of only 1 A/min was used, the patterns became indistinct due 
to a high background intensity, and no three-dimensional crystallites 
formed in the later stages of growth of the deposit. These effects were 
thought to be due to oxidation or attack of the lead by the residual 
oxygen and water vapour in the vacuum atmosphere. 
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3.2. Thallium 


Metallic thallium resembles lead in many of its physical properties, 
and it was found that the growth of a thallium film was closely analogous 
with that of the lead deposits described above. Thallium has a close- 
packed hexagonal structure and thin layers again gave rise to streak 
patterns with (001) thallium parallel to (111) silver and [100] thalium 
parallel to [110] silver, as shown in fig. 3 (Plate). It was found that the 
secondary diffraction pattern was much stronger than that which was 
observed from the lead films, but no satisfactory explanation for this is 
known. From the width of the primary streaks it was concluded that 
the deposit had formed islands of monolayer growth about 60 A in 
diameter. The closest distance of approach was measured relative to 
the silver and a value of a=3-31-0-02 A was obtained. This is about 4% 
lower than the value of a=3-45 A found for bulk thallium by x-rays. 
This contraction is of the same order as that found in the very thin lead 
deposits. 

As the thickness was increased beyond about 4 A, discrete spots became 
visible on the streaks, and rings due to hexagonal thallium appeared. 
When the thickness was increased up to about 100 A, a sharp pattern of 
slightly arced spots, together with sharp rings, was obtained. It was 
impossible to estimate the width of the oriented crystallites due to the 
arcing, which amounted to about +3° of disorientation. The width of 
the rings corresponded to a crystal size of about 100 A. 

The above results were obtained with a rate of deposition of about 
10 A/min, since at lower rates there was a deterioration of the diffraction 
patterns, presumably due to oxidation of the thallium. 


Shan LW: 


The results for tin were quite different from those described for lead 
and thallium. As the thickness of the deposit was increased, the pattern 
from the silver substrate weakened and two very diffuse haloes became 
visible, indicating that the deposit had an amorphous structure. The 
silver pattern gradually disappeared and then, quite suddenly, spots 
appeared along the centre diffraction line, as shown in fig. 4 (Plate), 
when the thickness of tin was about 7 to 8A. Deposits greater than 10 A 
in mean thickness gave rise to a sharp pattern of spots and rings as shown 
in fig. 5 (Plate), and the broad haloes, which can be seen in fig. 4 (Plate), 
completely disappeared. This sequence of patterns indicates that the 
amorphous layer of tin suddenly crystallizes at a critical thickness into 
both oriented and random crystallites. The sharpness of the spots and 
rings indicates a particle size of at least 200 A. 

The ring pattern was due to white tin which possesses a tetragonal 
structure with a=5-82 A and c=3-17 A and where the structure amplitude 
factor is given by | 


S={I+exp (wih) }[1-+exp {ai(hy+hg/2)}]. 
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An analysis of the spot patterns showed that there were three equivalent 
orientations of the tin, ie. (100) tin parallel to (111) silver and [010] 
tin parallel to any of the three equivalent (110) silver directions. The 
positions of the tin atoms relative to the substrate is shown 
diagrammatically in fig. 6. It is seen that there is only a small misfit 
along the [110] silver direction, since the distance between silver atoms 
is 2-89 4, which is almost half the distance between the tin atoms of 
5-82 A. In addition, there is an almost perfect angular fit since the 
[011] tin direction makes an angle of 28° 15’ with the [110] silver direction, 
and hence is almost coincident with the [2I1] silver direction. The 
relative misfit along the [001] tin direction is, however, —36%. 


Fig. 6 
Ag Cil2] 


Sn Cool] 
Ag [aii] 


Ag (110) 
Sn COI0) 


@ Silver atoms. ™ Tin atoms. 
The arrangement of tin atoms relative to the substrate silver atoms. 


There were also some spots on the centre diffraction line due to the 
orientation (110) tin parallel to (111) silver; it was not possible to 
determine whether there was any azimuthal orientation as there were 
no other associated strong spots. It is interesting to note that this was 
the orientation which predominated immediately after the transition 
from amorphous to crystalline tin, although the reason for this is not 
understood. 

3.4. Silver 

Deposits up to 50 A in average thickness produced no change from 
the original substrate diffraction pattern, ie. the orientation of the 
substrate was continued in the deposit, the surface remained smooth on 
an atomic scale (the diffraction spots still showed refraction), and the 
prominent Kikuchi pattern remained quite clear as shown in fig. 7 (Plate). 
As the thickness was increased beyond 50 4, some randomly oriented 
silver was formed and the Kikuchi pattern became much. fainter. 
Eventually the Kikuchi pattern disappeared completely and the propor- 
tion of randomly oriented silver increased until it predominated, the size 
of the crystallites being about 50 A. 
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3.5. Gold 


Thin layers up to about 8A in thickness were oriented in double 
positioning, it being quite easy to detect a marked change in the diffraction 
patterns, even though the pattern from the gold was superposed on that 
of the silver substrate. The refraction effects due to the flat surface of 
the silver disappeared and the diffraction spots became much broader. 
The width and height of the gold crystallites were estimated to be about 
40 A and 20 A respectively. 

When the thickness was increased beyond 8 A, broad rings due to 
randomly oriented gold appeared in the pattern and the Kikuchi pattern 
from the silver disappeared. Thicker layers gave rise to a ring pattern 
only, the sharpness of the rings increasing somewhat as the deposit 
built up. Measurements showed an increase in crystal size from 20 A 
to about 40 A. 


3.6. Nickel 


The experiments with nickel were not very satisfactory, since the 
molten bead of metal attacked the tungsten evaporator filament, causing 
it to burn out. Sufficient control over the rate of deposition was not 
therefore obtained, and the amount of tungsten evaporated with the — 
nickel may have been quite appreciable (see Heavens 1952). 

Oriented deposits showing double positioning were obtained however 
up to a mean thickness of about 2 to 34. The size of the crystallites 
was only about 25 A, and hence accurate measurements of the lattice 
constant were difficult to make because of the broadening of the spots. 
A value equal to that of bulk nickel was obtained but the error was at 
least +2%,. 

Thicker deposits showed only randomly oriented nickel with a small 
crystal size. 


3.7. Antimony 

Deposits of only a few angstroms in mean thickness had the amorphous 
structure similar to that found by Richter (1943), and Pinsker (1953, 
p. 244). The next stage in the growth of the deposit was noted when 
strong well-defined spots appeared on the centre diffraction line, as shown 
in fig. 8 (Plate). There were only very faint spots visible in other 
positions, and hence it must be concluded that there was a regular 
packing of layers of deposit atoms on the surface, but’ that there was 
little regular structure in each layer. As the thickness was increased 
still further, an extensive spot pattern appeared and the broad rings due 
to the amorphous antimony completely disappeared, as shown in 
fig. 9 (Plate). The spot pattern was independent of the azimuth setting 
apart from some very slight intensity changes, showing that the ~ 
crystallized antimony had a fibrous orientation with the pseudo-cube 
plane ((110) in the rhombohedral system) parallel to the silver (111) 
plane. This orientation is similar to that found by Geiling and Richter 
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(1949) for arsenic deposits on an amorphous base. The size of the 
antimony crystallites was estimated to be 254 high and 404 wide. 
Further experiments revealed that the fibre axis was not influenced by 
the angular position of the evaporator relative to the silver surface, 
as has been found in some cases (see Schulz 1949). 


3.8. Other Metals 


Of the metals examined, germanium and silicon formed only amorphous 
deposits, which is in accord with the work of Collins and Heavens (1952), 
and Hass (1948). Deposits of tellurium were completely amorphous in 
very thin layers but became partially crystalline as the thickness was 
increased (see Sakurai and Munesue 1952). Bismuth deposits were 
fibrously oriented with both the (110) and (211) planes (rhombohedral 
system), in contact with the silver (111) plane ; there was also randomly 
oriented material present. Deposits of aluminium, cobalt, chromium, 
iron, maganese and palladium were all randomly oriented, the particle 
size being about 20 Ain all cases. Zine and cadmium would not condense 
on the silver substrates even when high rates of evaporation were used. 
Similar effects have been found for these metals by other workers 
(see Cockcroft 1928 and Ehlers 1953). 


§ 4. Discussion 

The present experiments provide useful information about the initial 
stages in the formation of oriented overgrowths and indicate that either 
the type of growth mechanism proposed by Frenkel (1923) or the mono- 
layer mechanism may occur in any particular case. Thus, in the case of 
tin deposits, a transition from an amorphous to a crystalline state was 
observed when the average thickness of the deposit was about 10 A. 
This is consistent with the Frenkel model, in which it is assumed that 
deposited atoms still retain a certain amount of energy and move over 
the surface freely as a two-dimensional gas. As more atoms arrive on 
the surface, larger mobile aggregates form and at some stage spontaneous 
crystallization occurs. The problem is to predict where the crystallites 
will form and whether they will be oriented relative to the substrate. 
The crystallites may form at random points on the surface, or they may 
grow at the ends of dislocations or at gross faults in the substrate, as has 
been observed in the growth of copper on hot silver (see Kehoe et al, 1956). 
In the case of tin, both randomly oriented and oriented crystallites were 
formed, although their location was not known. 

The growth of antimony deposits was similar to that described above, 
although the antimony crystallites, once formed, had only a fibrous 
orientation with the plane of closest packing in contact with the silver 
surface. 

This mode of growth may also be described by the spontaneous 
crystallization of a supercooled liquid. Takagi (1954) found that thin 
layers of tin, lead and bismuth, 50 4 in thickness, showed remarkably 
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large supercooling effects; up to 192 degrees for tin, 157 degrees for lead, 
and 147 degrees for bismuth. He also showed that the degree of super- 
cooling increased with decreasing film thickness. Thus, in the case 
of tin which has a melting point of only 232°c, supercooling right down 
to room temperature might easily be obtained for the thinnest observable 
- deposits. 

Although there is some indirect evidence to suggest that this type of 
mechanism may occur to a limited extent in other cases, described below, 
it obviously does not apply generally to the same extent as in the growth 
of tin deposits. Thus, no amorphous or supercooled liquid state was 
observed in the initial stages of growth of bismuth, thallium or lead 
deposits even though these metals possess low melting points of 271°c, 
304°c and 327°c respectively. 

In the initial stages of growth of nickel, gold and copper (see Newman 
and Pashley 1955), oriented crystallites were formed. It was found that 
individual crystallites did not continue to grow as more material was 
deposited, but that the density of crystallites increased. Similar results 
have also been found by other workers using different materials (see 
Schulz 1952). These observations are in accord with a spontaneous 
nucleation process, although the average thickness when crystallization 
occurs would have to be very small, since no amorphous deposits are 
observed. The reason why the crystallites do not continue to grow is 
not known. In the case of the copper deposits, it may be argued that 
a thin surface film of oxide forms on the crystallites and prevents their 
further growth, but this argument cannot be applied to gold deposits. 
These observations explain in a simple way why deposits of these materials 
rapidly become randomly oriented when the thickness is increased beyond 
a few angstroms. When the density of crystallites becomes such that 
an appreciable proportion of the surface is covered, growth must occur 
on top of the initial crystallites. In general, there will be regions of 
misfit between neighbouring crystallites, since the lattice constant of the 
deposit bears no simple relationship to that of the substrate, and the 
crystallites grow at random points on the surface. Random growth 
would be expected from these regions of misfit. These results emphasize 
the need to examine the thinnest possible films when testing whether or 
not epitaxy occurs in a particular case. 

Another type of growth occurred in the case of lead and thallium 
deposits which grew in the form of monolayers. Three-dimensional 
crystallites were only observed when the average thickness of the deposit 
was increased beyond about 4A. Thicker deposits again showed a 
mixture of oriented and randomly oriented crystallites. The randomly 
oriented material may have formed by the mechanism given above. 
The growth of silver on itself must also have proceeded by a monolayer 
mechanism because of the observed refraction effects from layers up 
to 60 A in thickness. The growth of silver was different from that of the 
other face-centred cubic metals, since it grew in a strictly parallel 
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orientation, in contrast to the orientation in double positioning in the 
other cases. 

The reason for the non-orientation of certain metals such as cobalt, 
iron and palladium, may be due to low mobilities of atoms of these 
materials on the silver surface. This would be substantiated by the 
small crystal size observed for these deposits. Oriented growths of iron 
and palladium on a copper single crystal have, however, recently been 
found by Haase (1956). 

The parameter of relative misfit between deposit and substrate is 
often regarded as an important factor in the formation of oriented 
overgrowths and this factor has been found to influence the nature of 
deposits in some cases. Thus, Schulz (1951, 1952) found that alkali 
halides grew as monolayers for small values of misfit with the substrate 
(less than 10°, on mica, and less than 20°, on other alkali halides), 
whereas nucleated growth occurred for larger misfits. Such an effect 
has not been observed in the present work. On the contrary, gold, 
nickel and copper grew as nucleated deposits, although the misfit is 
less than 14° in each case, while lead and thallium grew in the form of 
monolayers where the misfits are 21° and 19° respectively. 

The lattice spacing of the deposit in the initial stages of growth is also 
of importance as a result of the theory of Frank and Van der Merwe 
(1949 a, b) and previously reported cases of pseudomorphism (for a review 
see Pashley 1956). If the deposit grows in the form of oriented three- 
dimensional crystallites which increase in density as the total amount 
of material is increased, then the growth of an individual crystallite 
cannot be followed. Thus, the lattice spacing of the atoms in the plane 
in contact with the substrate cannot be measured. Pseudomorphism 
could only be detected if the overgrowth builds up atom by atom in the 
form of a monolayer. Lead and thallium have grown in this way in the 
present experiments. The results show that the deposits had their own 
characteristic lattice constants in the earliest stages of their growth. 
The spacings found were somewhat modified from those of the bulk 
metals, but this may be expected as a result of the two-dimensional 
character of a monolayer. There was no evidence to support the theory 
of Frank and Van der Merwe (1949 a, b) which predicts that oriented 
monolayers are constrained to fit the substrate lattice. These two cases 
probably afford the best evidence on monolayer growth, although Schulz 
(1951, 1952) has also observed monolayer growth with alkali halides. 
The metal deposits seem to be somewhat superior than those due to 
the larger degree of misfit which results in a better resolution in the 
diffraction patterns. Also, no difficulty is experienced due to charging 
effects which cause a blurring of the patterns. 

Finally, it should be remembered that the vacuum conditions under 
which such work as that described above was carried out, are far from 
satisfactory. It is well known that surfaces are covered with a mono- 
layer of gas when the vacuum is only 10° mm Hg. Ideally, experiments 
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should be performed in a vacuum of at least 10-® mm Hg when the rate 
of bombardment of the surface by the gas molecules would be less than 
the rate at which the deposit atoms arrive (1 A/min). Results obtained 
under these conditions may be vastly different from those described above, 
with regard to the size of crystallites and the degree of orientation. 
Unfortunately a vacuum of this order is not at present attainable in 
a standard metal apparatus. 
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ABSTRACT 


The interaction between a positive ion and a neutral atom is calculated 
from a statistical model which treats the electrons as a Fermi gas and 
considers the electrostatic forces between it and the nuclei, as well as the 
kinetic and exchange energies of the gas itself. This is applied to the case 
of Ca*~* ions and neon atoms and of K* ions and argon atoms. Theformer 
is extended to consider also neutral calcium atoms and neon atoms. 

The elastic scattering of these ions is calculated by the classical and 
also by the quantum method over an energy range up to about 15 kev 
using these interactions. The results have been applied to problems 
of meteor trails (Massey and Sida 1955) and the spread of a beam of 
ions due to multiple scattering. 


§ 1. INTRODUCTION 


ALTHOUGH the scattering of electrons and protons has been explored 
theoretically, little has been accomplished for heavy atoms. These are 
important, for many processes, both in the laboratory and in astrophysics, 
involve the scattering between heavy atoms and ions. For example a 
knowledge of the behaviour of a beam of positive ions in a neutral gas 
is often required and such questions as its spread depend on the amount 
of scattering. In astrophysics the processes of excitation and ionization 
of meteor atoms are brought about by collisions with air molecules and 
atoms and using the methods of this paper these have been discussed by 
Massey and Sida (1955). The first requirement of the scattering problem 
is to obtain an accurate representation of the inter-atomic forces between 
the particles. 

The Hellmann—Feynman theorem (Feynman 1939) indicates that the 
forces on the nuclei of a molecule as determined from quantum mechanical 
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potential energy surfaces are essentially those arising from classical electro- 
statics and a knowledge of the electron probability density. Thus the 
inter-atomic forces will depend on obtaining this electron charge density 
and determining the electrostatic forces which arise from the nuclear 
charges and the inter-penetration of the electron shells. Most methods 
which have been used have been based on this theorem but the difficulty 
is to represent the electron density and the part played by the Pauli 
Principle. To overcome these problems it has been usual to postulate 
a form for the repulsive part of the potential as an exponential or a 
polynomial in inverse powers, and to obtain values for the parameters 
thereby introduced from considerations of such properties as viscosity, 
the virial coefficients and crystal structures (Massey and Burhop 1953, 
chap. VII). 

In this paper we do not postulate any particular form for the interaction 
between two atomic systems but make certain assumptions by which 
the electrons of the atoms concerned are treated statistically. The 
statistical model employed here is substantially that used for solid 
state investigations by Lenz and Jensen (1932) and refined and widely 
applied by Wigner and Seitz (1934) and by Gombas (1949). The electrons 
are treated in each particle as a free Fermi gas and we consider the change 
in the total energy of the two systems caused by overlapping of the two 
gases and also the electrostatic forces arising between the nuclei and 
these clouds.. Thus one expects the method to be most applicable to 
particles with (a) many electrons and (b) closed shells (e.g. Catt, Kt 
A, Ne). There is a further restriction which may be important at high 
energies in that the particles must not approach each other too closely 
as no allowance is made for distortion of the original cloud distributions. 

At large distances of nuclear separation the chemical repulsive forces 
are replaced by the attractive Van der Waals forces. To a first approxi- 
mation these can be expressed as a single inverse power but the present 
method fails to give these terms, as indeed most other methods have 
failed. The attractive forces however give no contribution for the 
scattering of heavy particles above about 50 ev and hence are neglected 
in this paper. 

In §2 the interaction expression is derived and in §3 applied to the 
two cases of calcium ions with neon atoms, and potassium ions with 
argon atoms. In §4 these results are applied to calculations of elastic 
scattering and finally the problem of multiple scattering of a positive ion 
beam is also discussed. 


§ 2. THE INTERACTION BETWEEN ATOMIC SysTEMS 


The outline of the method used in calculating the interaction between 
the systems is given in the paper by Massey and Sida (1955). The first 
formula is eqn. (8) of that paper which can be written in the form 


Vrar [Wat Wee te W gt get gene, em 
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W,=0[Z,2,—{Z2,N.4+Z,N}+N,Ng], 
W = e'[Z,N o90(7) + ZN 19,(7r)], 
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Ir—ral 
V(r) is the interaction energy between the two systems when at a distance 
r apart. Z,, Z, are the nuclear charges and N,, N, the numbers of 
electrons in the respective systems. 4, and 4, are the electrostatic 
potentials due to the electron clouds in each atom and 


d{t)==— a t—g(r)}. 

A and B are the two constants 
A=0-5312a,e"3, B=0-3 17 764% 

The first term W,, is the covlomb interaction of the positive and negative 
charges and is zero if either particle is neutral. W, arises from one 
nucleus of charge Z,e being penetrated by the charge cloud of the second. 
The cloud no longer screens its nucleus completely and the penetrated 
nucleus is thus repelled by the effective nuclear charge, V.g,(r), of the 
approaching atom. W, arises from those parts of the charge clouds 
which overlap and hence no longer repel each other : the unsymmetrical 
form is due to simplifying by Green’s theorem. 

The interaction V(r) can be evaluated if suitable forms for Ng(r) and 
V?d(r) are obtained for each atom. These represent the effective nuclear 
charge and charge density respectively and can be obtained from either 
the Fermi-Thomas—Dirac model or the Hartree-Fock field. The latter 
was chosen as being more accurate since the results are available for the 
atoms investigated. The two cases investigated were Ca** ions with 
neon atoms and K* ions with argon atoms. The integrations W,, W, 
and W, were performed numerically and were rather laborious. The 
values of Ng(r) and Vd(r) for each atom were taken from the Hartree— 
Fock calculations for calcium (Hartree and Hartree 1938), potassium 
and argon (Hartree and Hartree 1935) and for neon (Brown 1933). In 
terms of Hartree’s results 


where Z, and U(r) are tabulated. 
We also investigated the effect of the two valence electrons, (*S)?, 
of calcium and V(r) was calculated by amending the values of Ng(r) 


and V?4(r) accordingly. 
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§ 3. RESULTS OF THE CALCULATIONS 


From the form of V(r), eqn. (1), being the sum of five terms it follows 
that there will be a certain amount of cancellation and this coupled 
with the tedious numerical calculations is rather a disadvantage to the 
method. Tables 1 and 2 give the results including values of the respective 
terms. 

Table 1. Interaction of Ca** Ions with Neon Atoms 


r +W, —W, | +W, | =—W,.| 2W, | V=Qp2R, 
0-5 | 62:20 | 36-89 | 23-06 | 1-80 46-57 93-14 
0-8 | 30-40 | 23-98 | 21-03 | 2-05 25-40 31-75 
1-0 | 18-60 | 18-32 19-23 | 1-94 17-57 17-57 
1-4 7-20 | 10-12 12:39 | 1-60 7-87 5-62 
ee a 3-51 5-71 | 0-94 3-06 1-53 

(3-66) | (6-11) | (1-09) (3-16) (1-58) 
an 9.299 | 0-486 | 0-947 | 0-257 0-424 0-141 

(0-762) | (1-578) |(0-422) | (0-613) (0-204) 

; 0-0635 | 0-121 | 0-055 0-0305 0-0076 

4-0 0-032 

(0-296) | (0-711) |(0-228) | (0-219) (0:0548) 
Le 0-012 | 9:0285 | 0-043 | 0-025 0-007 0-0015 

(0-207) | (0-552) |(0-180) | (0-178) (0-0395) 


Figures in parenthesis refer to amendment due to (48)? valence electrons for 
calcium. Unit of r is a) and of V(r) it is e?/ay. 


Table 2. Interaction of K* Ions with Argon Atoms 


r +W, —Ws, We 2 We =W; |V=(n2W, 
0-5 | 102-12 63-99 | 38:30 | 2-52 | 73-91 147-82 

0-8 | 65-550 | 48-756 | 24-148 | 0.2300 | 38-552 48-19 
1.0 | 47-690 | 40-530 | 23-897 | 9-432 | 28-626 28-625 
2-0 7-077, | 11-079,-|°10:901 | 1-710 | 1890 2-594, 
3.0 0-9025 | 2-0120 | 3-1699 | 0-5725 | 1-308 0-436 
4-0 0-1045 | 0-3040 | 0-7100 | 0-2327 | 09-2778 0-0694, 


5-0 0-0125, 0-0400 0:1389 | 0-0640 0-0474, 0-0095 


a 


Unit of r is a and of V(r) it is e%/ay. 
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The exchange term W has little effect for <1 but is of considerable 
importance at large distances. It can be shown that as r—> 0, W,and W, 
each tend to constant values. 

For values of r>4, the interactions all tend to zero exponentially 
so that no attractive force is given. This forms a limit to the validity 
of the method and it would be unwise to extrapolate much further. 
So that in all applications we must ensure that any contribution to 
scattering from parts of the potential curve for r>4 is quite small ; this 
forms a lower limit to the energy and angle at which scattering can 
be investigated. It is impossible to carry out calculations of V(7) at much 
larger distances since there are insufficient figures given in the Hartree- 
Fock results to allow for the cancellation which always occurs. 

Actually Gombas (1943) introduced a further term (Wg), similar 
to W, to take account of electrons with anti-parallel spin, a polarization 
or correlation term, which he evaluated for the Fermi~Thomas-—Dirac 
field. This he estimated increases the electron density in the marginal 
region of the atom by about 45% in his model but the exchange term 
by only 13%. Wigner (1934) also introduced a correlation term in his 
solution to the Hartree-Fock field with considerable improvement to his 
results which depended on long range interactions. However the work 
involved in such an amendment to the present scheme does not seem 
justified. 


$4. ExLastic SCATTERING AND Momentum Loss Cross SECTION 
The classical expression for the differential elastic scattering cross 
section should be quite valid for the heavy mass particles with which 
we are concerned down to very small angles and low energy. This is 
given by 


I(@) sin 6 dé df6=p z sin 6 dé dd 


dp 
dO 


where 


— 
bo 
ae 


iv.) 

pn —2p | dr/r{r2(1— V/E')—p?}*?, 
To 

p is the impact parameter, 7) is the largest singularity of the integrand. 

(Massey and Burhop 1953, p. 373). From this expression the momentum 

loss cross section is given by 


Q (E)=27 | a eos OVI) sin. Odd: © wr S. -w= (8) 


This converges since the (1— cos 6) factor considerably decreases the 
effect of small angle scattering which according to (2) increases to 
infinity. Results for collisions of K* ions with energies ranging from 
100 ev to 30000 ev with argon atoms at rest are given in table 3 (Q5) 
and in figs. 1 and 2 /(@). The only experimental observations available 
for comparison are those of Rouse (1937) which are rather limited in 
angular range and are not absolute. As far as it can be checked there is 


good general agreement. 
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Cross sections for collisions of Cat + ions with neon atoms have already 
been published by Massey and Sida (1955). 

The quantum mechanical method gives the total elastic cross section 
and momentum loss cross section in terms of the phase shifts as follows 


(Mott and Massey 1949). 
4 A 
Q(E)= 5 D(2n-+1) sin? n, peter sereamciie re (C8 


Q(B) =o Sim) sin? (Maa) ss + 8) 


The B.W.K. approximation 


co ¢ T\2>) 1/2 oo 1)277 1/2 


where 7, and 7p are the roots of the integrands, was used for large values 
of 7,, and the Born approximation 


a 2m 


ne — gar) VOW ernlbrPrdr - (7) 


for small values up to a radian. It was found that for these heavy 
atoms the Born approximation was quite reliable up to this value. 


Table 3. Momentum Loss Cross Section for Collisions of K* Ions 
with Argon Atoms 


Energy Ga Energy Qa 
(ev) (units 7a”) (ev) (units 79") 
180 4-21, 4000 0-66 
360 3-05; 10000 0-300 
1000 1-68 20000 0-157 
2000 1-08 30000 0-107 


For calcium ions at 200 ev and potassium ions at 30 kev the momentum 
loss was calculated by both methods with good agreement in both cases. 

There was no difference in the calculations over a range 100-2000 ev. 
for momentum loss or large angle scattering of Cat*+ ions and neutral 
calcium atoms by neon atoms; the difference was only in small angle 
scattering (<20°) and is shown in the total elastic cross section calculated 
from (4) which was 377a,? for Ca** but 1027a,? for the neutral atoms. 

The difficulty of experimental verification is indicated in the results 
for K* ions at 30 kev for which the quantum method gives 5x 109a,” 
for zero angle scattering /(0) while the classical method at 0-005 radians 
gives for [(@) 8x 104a,?; this is a difference of 5 orders of magnitude 
in about 20 minutes of arc, for with such heavy particles classical scattering 
should be valid down to about 10 minutes of are. 
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For large energies, above 10 kev, the results for ions indicate that the 


total momentum loss cross section varies approximately inversely with 
the energy. 


Fig. 1 


(1 (8) in units a2 per steradian) 


1(8) sin 8 


——_ 


' A 2 . = 
Angle of scattering (radians) (c.ofm. system) 


Elastic scattering of K+ ions by argon atoms. 
Energy 2-30 kev (lab. system). 


Fig. 2 


1 (8) sin® 
(I(Q)in units a,2 Per steradian) 


! 2 3 
Angle of scattering Cradians)(c. of m. system) 


Elastic scattering of K+ ions by argon atoms, 
Energy 180 ev—2000 ev (lab. system). 


§ 5. SPREAD OF AN Ion BEAM DUE TO MULTIPLE SCATTERING 


A narrow beam of positive ions will be spread by multiple scattering 
with gas atoms; the amount of spread may be calculated as a function 
of the path length and pressure of the residual gas through which the 
beam passes. This is particularly important for isotope separation 
and similar experiments. From the above work it may be very simply 
evaluated to a good approximation. . 

Consider a particle incident in the z-direction on a scattering material. 
We suppose that in its collisions with the scatterer the particle loses 


3F2 
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a negligible fraction of its incident energy. Then, according to Fermi, 

the probability that, after penetrating a distance z into the scatterer, 

the transverse displacement of the particle is x and projection of its 

direction of motion on the wz plane makes an angle 7 with the z-axis, is 
4/3 x etme eee a ee 

CE Re balled — -(7?—3n-+ > eS 

Wes 52) dm dij 5 CXR | 9 eg) a 


The mean free path is given in terms of a cross section @, by 


ASIN 2s a 
where N is the number density of particles and 
rot hot 
Qa= | | (1— cos #)1(0) sin 0 d0dé . . . . (10) 
0/0 


where /(@) sin 6 dé d¢ is the differential cross section for elastic scattering 
of the particles by the scatterer. This assumes that the scattering centres 
are infinitely massive. In our case in which the incident ions and scattering 
atoms are of nearly equal mass we must replace (10) by 


2m 
Qd= |_| (1— cos 40)1(9) sin0d0dé . . . . (1) 
070 


which we shall call the effective diffusion cross section. 

Scott (1949) expresses the formula (8) in terms of the angular 
displacement % and the direction ¢ to the chord of the position of the 
particle. Thus fa/z and d=n—w. Hence in eqn. (8) 

4/3 X AG : 
Wohi bs ab dp YP exp{— 2 2 db-+8)} adh ds.. (12) 


ee 


Table 4. Effective Diffusion Cross Sections Q ,° for K* Ions in Argon 


Energy (lab. system, kev) Q (749°) 


30 0-030 
20 0-045 
10 0-086 

t 0-193 


We transform to a dimensionless parameter 
p=\e—1NeO eee 2 ae ls) 
which is the ratio of the mean free path to the total ‘path traversed. 


Then on integrating over the angle ¢ we get the total probability of an 
angular displacement in one plane between y and 4-+dib to be Wb; 2) deb 


where 3y\ 1/2 ; 
Pp 3 
Wb; z)= (3) exp Graak - 7 ee ley 


For potassium ions in argon the effective diffusion cross section 05° 
is found to have the values given in table 4. 
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Calculations were then made to obtain the distribution function (14). 
in a number of different cases of potassium ion beams in argon to: 
investigate the effect of varying the total path length, the pressure, and 


Hig. 3 


W(y, z) 


0.04 0.02 ° 0.02 0.04 
Deviation ¥ in radians 
Multiple scattering of K+ ions in argon. Variation with path length (z)- 
Curve 1 z=20cm. Curve I] z=40cm. Curve III z=100 cm. 


(Pressure of gas 1-3 10-° mm. Cross section Q ;’=0-117a,?.) 
Incident Energy 8 kev (lab.). 


Fig. 4 


Wy, z) 


0.04. 0.02 re) 0.02 0.04 
Deviation ¥ in radians 


Multiple scattering of K+ ions in argon. Variation with pressure (p). 
Curve I p=1-:3X10- mm. Curve II p=1-3 x 10-° mm. 
Curve III p=1:3 x 10-4 mm. 


(Path length z=20 cm. Cross section Q;=0-117a,?.) 
Incident energy 8 kev (lab.). 


the cross section (i.e. incident energy). The results are shown in figs. 3, 4 
and 5. The considerable broadening of the beam for long path lengths, 
high pressures and low energies is clearly shown and indicates the degree: 
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to which these factors have to be taken into account in retaining a narrow 
beam, or in the case of isotopes in obtaining a high yield of separation. 
It must be remembered that no allowance has been made for the energy 
loss in small angle elastic scattering. This will also increase the 
broadening particularly in isotope separation. Its importance will be 
examined in a later paper. 


Fig. 5 
os 1 
= 
= 
> 
0.02 O 0.02 


Deviation ¥ in radians 
Multiple scattering of K+ ions in argon. Variation with cross section (QQ). 


Curve I Q,=0-0557a,?, energy 16 kev (lab.). Curve Il Q,=0-ll7a,”, energy 
8 kev (lab.). Curve III Q,=0-16za,", energy 5 kev (lab.). 


(Path length z=20 cm. Pressure p=1-3 x 10-° mm.) 


§ 6. EXTENSION TO VERY Low ENERGIES 


The present method has been shown to be suitable where the scattering 
is dependent only on the repulsive forces between the particles. For 
very low energies of several hundred degrees absolute, the attractive 
mteractions associated with large inter-nuclear separations are of prime 
importance. Therefore for problems of this nature a polarization energy 
inversely proportional to the fourth power of this separation, can be 
postulated. It is not altogether clear under what circumstances the 
scattering is independent of the repulsive interaction. A calculation 
of the mobility of potassium ions in argon atoms at 293°K gave a value 
20% below the experimental result when only a polarization potential 
was used. The result differed even further from the Langevin theory and 
furthermore the variation with temperature was also very different. 
This suggests that the region of the interaction where it reaches a minimum 
and where it changes from positive to negative values must be taken into 
account for low energy scattering. At present there is no obvious way 
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of fitting a polarization potential to the statistical results and it would 
seem that this is required so that we have a continuous curve. The 
most likely method would be to investigate the potential forms required 
to give the experimental mobility results in the same way that the viscosity 
has been used for calculations on the interactions between neutral atoms. 
The wide range and accuracy of the experimental data make this possible 
and such a programme is obviously desirable from a theoretical standpoint. 
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ABSTRACT 


The plastic deformation of single crystal bars of silver chloride in 
tension and bending has been studied by using the photoelastic effect. 
It is found that in tension the plastic deformation is not uniform over a 
given cross section but is greatest near the surfaces. In bending the 
plastic deformation is found to propagate longitudinally along the bars 
rather than transversely. The distribution of birefringence, and hence 
that of stress, has been measured and the result interpreted by means of 
a simple dislocation model. It is concluded that, in silver chloride, the 
ratio of the stress necessary to propagate dislocations to the stress 
necessary to generate them is less than about 0-6. The shear stress on 
the glide plane necessary to propagate dislocations is approximately 
5-2 x 10° dynes/cm?. 


§ 1. INTRODUCTION 


WHEN a crystal is not fully annealed the dislocations within it produce 
long range elastic strains. In a transparent cubic crystal such strains 
are visible as distributions of birefringence produced by the photo- 
elastic effect. Detailed photoelastic observations of silver chloride 
under a low power microscope (Nye 1949 a, b) have shown that, after 
plastic deformation by glide, the internal stresses are of two kinds (other 
cubic crystals show similar effects). (A) Each glide packet is subject to 
pure elastic bending; this gives a ‘ saw-tooth’ distribution of stresses 
with a period equal to the spacing of the glide zones. (B) Superimposed 
on the saw-tooth distribution is a more slowly varying stress ; it is this 
second, larger scale, kind of stress which is the subject of the present note. 


§ 2. EXPERIMENTS 
2.1. Plastic Hutension 
Single crystal bars of silver chloride have been prepared with rectan- 


gular cross sections of dimensions 1x1 to 1x5 mm and lengths from 
30 to 45mm. The final treatment of the surfaces was a polish produced 
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by agitation in a strong solution of ‘hypo’ for a few minutes. The 
orientations were variable. When a fully annealed crystal is plastically 
extended and the load is removed a photoelastic pattern is seen between 
crossed nicols which corresponds to stresses of type (B). It always has 
the general form shown in fig. 1(a)+. In this figure two neutral lines are 
visible parallel to the length of the specimen. (The path differences in 
this, as in nearly all our experiments, are less than a quarter of a wave- 
length and the extinction angle is nearly the same at all points of any one 
specimen.) The neutral lines also show up well in a circular polariscope. 
The variation of the birefringence across the specimen, which gives rise to 
the photoelastic pattern, is shown by the curve in fig. 1(b). Since we may 
assume that the birefringence is proportional to the stress, the curve may 
equally well be taken to represent the distribution of the (longitudinal) 
internal stress o in the specimen. A photoelastic pattern of the same 
general type is seen when the bar is viewed through the other pair of 
opposite faces. The sign of the effect is such that the outer parts of the 
bar are in compression, while the inner core is in tension. It seems, 
therefore, that during the plastic extension the outer layers have deformed 
plastically more than the inner core ; for when the load is removed the 
outer part would then be thrown into compression, as is observed. One 
is thus led to conclude tentatively that more dislocation sources are 
active at or near the surface than in the interior. 


2.2. Plastic Bending 


When a similar specimen is bent elastically by applying equal and 
opposite couples to the ends, so that the upper side is convex, the photo- 
elastic pattern seen between crossed nicols is as shown in fig. 2(b). The 
corresponding distribution of birefringence, and therefore, we suppose, of 
longitudinal stress, is drawn in fig. 2(a). The bar thus shows the normal 
linear variation of stress to be expected within the elastic range, and if 
the couples are removed the birefringence disappears. (As will be 
explained later, the specimens from which the photographs of fig. 2(b, c, d) 
are taken were actually partly in the elastic and partly in the plastic 
range ; the section shown in fig. 2(b), however, was completely elastic and 
serves well enough for illustration.) When the bar is bent into the plastic 
range the photoelastic pattern changes to that shown in fig. 2(d) and the 
corresponding distribution of birefringence and stress is then as shown in 
fig. 2(e). If now the couples are removed from the ends the photoelastic 
pattern changes to fig. 3(a) and the corresponding stress distribution to 
fig. 3(b). There are now three neutral planes and the resultant distri- 
bution gives no resultant force or couple (bending moment) over a Cross 
section. Figure 3(b) is the same as fig. 2(¢) except for the addition of a 
linear elastic stress distribution ; we therefore conclude that there is no 
significant reverse plastic deformation when the external couples are 


removed. 
+All figures are plates. 
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The sign of the residual stress in these bending experiments is again con- 
sistent with the hypothesis that more plastic deformation has taken place 
near the surface than in the interior—but this is hardly surprising, since 
the outer layers were the most highly stressed during the bending. What 
is surprising is (i) the very high stress gradient across the neutral plane in 
fig. 2(e), as shown by the sharpness of the central dark line in fig. 2(@), and 
(ii) the uniformity of the stress in all other regions, except possibly in the 
regions very close to the edges. The transition from a uniform tension 
co, to a uniform compression —o,,, as judged by the change in sign of the 
birefringence, occurs over a distance which is typically about 0-1 mm, 
The actual transition may be even sharper, because it is probably blurred 
by overlapping of successive layers of crystal normal to the direction of 
viewing. The birefringence is uniform outside the transition region to 
within 5%. 

One naturally asks what are the intermediate stages by which the 
pattern of fig. 2(d) is produced. The conventional view of plastic bending 
is that the deformation starts at the more highly stressed outer surfaces, 
and then two plastic regions spread inwards towards the neutral plane, 
one from each side of the bar. But in our specimens this was not the case. 
Instead of propagating transversely, the deformation propagated longi- 
tudinally down the bar. It began at one end near the knife-edges used 
to produce the couples, and progressed down the bar at a rate which 
could be controlled by the rate of bending. Figures 2(b), (c) and (d) show 
three sections of bars in which the plastic zone had reached an inter- 
mediate position and with the external couples still applied. The section 
in fig. 2(b) is still in the elastic range, as already described. The section 
in fig. 2(d) is plastically bent. Figure 2(c) shows the transition region ; 
its length is rather less than the width of the specimen. 

To find numerical values for the stresses we need to know the photo- 
elastic constant, but since cubic crystals are not photoelastically isotropic 
the appropriate photoelastic constant depends on the crystal orientation. 
We have therefore measured the constant for each of our specimens by 
applying known bending couples in the elastic range and measuring the 
birefringence at known distances from the neutral plane. This enables 
us to deduce the longitudinal stress at any point in a plastic or partially 
plastic bar. In a partially plastic bar equilibrium considerations show 
that the maximum stress o,, in the elastic part must be 3 times the 
uniform stress oc, in the plastic part. This is verified by our experiments. 
Values of o,, in five different specimens were 12-7, 17-6, 19-4, 11-7, 
18-8 10° dynes/em?. Since deformation appeared to be by glide on 
a single set of planes we have calculated the corresponding shear stresses 
on the glide plane (not resolved along the glide direction, since this was 
not known). The results are 5-7, 8-5, 9-7, 5:8, 9-0x 106 dynes/cm2 ; 
mean=7-8 x 10° dynes/em?. The shear stress on the glide plane 
in the plastic region is } of these values and thus has a mean of 
5:2 10° dynes/cm?. 
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§ 3. Tenrattve Mopen 


Our tentative explanation of the behaviour of the bars in bending is 
briefly as follows. We may first try to explain how the distribution of 
stress shown in fig. 3(b) can be maintained. As a start, consider a row of 
negative edge dislocations arranged along the neutral plane with Burgers 
vectors parallel to the neutral plane and with their lengths parallel to the 
axis of bending. They would produce-a stress ‘ discontinuity ’ on the 
neutral plane of the same sign as that observed, that is, tension above 
and compression below (Nye 1949 b). But since the extra half-planes 
would be in the lower half of the crystal the bar would be bent with the 
upper surface concave ; this is the opposite curvature to that observed. 
Such a model therefore fails. 

We can, however, produce both a curvature and a stress discontinuity 
of the correct sign by taking the following related model. The crystal 
is filled with a macroscopically continuous distribution of positive dis- 
locations. These produce a curvature of the correct sign, but if arranged 
suitably (Nye 1953) they give no macroscopic internal stress. We now 
remove those (positive) dislocations which are situated within a narrow 
band containing the neutral plane. This is equivalent to adding a 
narrow strip of negative dislocations and produces a stress discontinuity 
of the correct sign. According to this model, then, the observed stress 
distribution is produced not so much by the presence of dislocations as 
by an absence of them. 

In the real crystal bars the Burgers vectors of the dislocations 
responsible for the bending are not parallel to the neutral plane. We 
deal with this by superposing on our model a distribution of edge dis- 
locations with Burgers vectors perpendicular to the axis of the bar. This 
second distribution produces no long-range stresses but accomplishes the 
necessary adjustment of the resultant Burgers vectors. There will also 
doubtless be present in the real crystal bars other dislocations apart from 
those responsible for the bending, but since their arrangement is such as 
to produce no long range stresses they make no contribution to the photo- 
elastic effect. 

An explanation for the longitudinal propagation shown in fig. 2 may 
be derived from a recent paper by Read (1957). Read has calculated the 
bending moment versus curvature curve for plastic bending by using a 
simple dislocation model. In his model o, is the stress necessary to 
generate dislocations and oy is the stress necessary to move them through 
the crystal. The curve obtained depends on the ratio o9/o, and is such 
that the bending moment goes through a maximum for all values of o/c, 
between 0 and 1. However, when o,/c, is greater than about 0-67 the 
maximum is very flat, and would easily be destroyed if a little strain- 
hardening were introduced. On the other hand, when o)/c, is less than 
0-6 the bending moment goes through a sharp maximum at the elastic 
limit. This, from our point of view, is the crucial fact, for we should then 
expect instability in a bending test ; when deformation started at a 
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certain place it would continue there preferentially. (One is reminded 
of the behaviour of mild steel in tension.) The fact that the deformation 
is non-uniform along the length of the silver chloride bars can therefore 
be interpreted to mean that in this crystal oy/c, is less than about 0-6. 

Apart from the instability just mentioned, which is a logical consequence 
of Read’s model, our results for the distribution of stress in the plastic 
region shown in fig. 2(e) are in general agreement with Read’s theoretical. 
expectations. In particular, the uniform stress c, that we observe may 
be identified with the theoretical stress oy. (For completeness it should 
also be said that according to Read’s model, the stress in the regions very 
close to the upper and lower surfaces in fig. 2(e¢) would not be o, but would 
rise linearly from oy to the value co, on the surface. Unfortunately, how- 
ever, the viewing conditions in these narrow regions very close to the 
surfaces are not good enough to enable us to say whether or not they 
behave as the theory predicts.) 

We have seen that the stress c,,, and therefore the stress o) necessary to- 
move dislocations through the crystal, is 5-2 10° dynes/em? when 
expressed as a shear stress on the glide plane. Our experiments give no 
precise information as yet about the value of the shear stress necessary 
to generate the dislocations beyond the fact that it must be greater than 
about 1/0-6 times this value, that is, greater than 8-7 x 10° dynes/cm?. 

More detailed measurements are being made of the phenomena both 
in pure tension and in bending, and a full account will be published. 
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In all previous investigations of the absorption of slow negative 7-megons 
by the constituent elements in nuclear emulsions, the separation of the 
observed events into those involving the nuclei of the light elements 
(carbon, nitrogen and oxygen), and of the heavy elements (silver and 
bromine) has been achieved largely by means of the ‘ «-particle method ’. 
In order to effect a separation of the events in light and heavy nuclei 
independently of the «-particle, and other indirect methods, an experi- 
ment has been performed in which 7-mesons were allowed to pass into 
emulsions of different compositions. 

Uford normal G5 and dilute G54 emulsions were exposed simul- 
taneously and symmetrically to z-mesons obtained by means of the 
Glasgow 300 Mey synchrotron. The stars induced by the 7~-mesons 
which stopped in the emulsions were subdivided according to the number 
n, of secondary prongs. 

The numbers of stars with n >1 have been corrected for the efficiencies of 
the observers for different types of event in normal and dilute emulsions. 
The relative numbers in dilute and normal emulsions have also been 
normalized to allow for differences in the areas scanned in each type of 
emulsion, for the differences in stopping power and for slight differences 
in thickness. The final normalization factor is in agreement with that 
obtained by comparing the numbers of positive 7-mesons stopping in the 
two types of emulsion. The prong distributions were obtained with the 
conventions (a) that only tracks of length +5 were counted as prongs 
and (6) that all tracks having a definite direction were counted as prongs. 
The results presented are those obtained using convention (5). 

If o, and o;, are the cross sections for the production of a particular 
type of event in heavy and light elements respectively, then the equivalent 
numbers of such an event in the normal and dilute emulsions are given by 

Mrormat C 2°060y-+ 2-870, 

Nanute © 0°920 q+ 4320, 
where the numerical coefficients are proportional to the numbers of 
heavy and light nuclei in the two types of emulsion. These equations 
may be solved to give oj and o, and hence the percentages of the various 
types of events in heavy and light elements for the normal emulsion. 

In order to decrease the statistical errors in oy, and o,, for the normal 
emulsion we have added to the data of the present experiment for this 
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type of emulsion those of all the previous workers who have used a 
criterion similar to ours for the minimum prong length. Thus well 
determined values are obtained for Nora: 20d the statistical errors in 
the final results arise from the more limited numbers of events, from the 
present experiment alone, which contribute to N ainute- 

The distribution in prong number for the heavy elements is shown in 
fig. 1(a) as a percentage of all stars in heavy nuclei, while fig. 1(6) shows 
the corresponding results for light nuclei. The errors shown are standard 
deviations and include the errors in the efficiency determinations. In 
order to compare the results with those obtained by other workers the 
numbers are presented in the table as percentages of the total number of 
7 stars in nuclei other than hydrogen in a normal emulsion. The 
results of the present experiment offer direct confirmation of the results 
obtained by the more indirect methods. 


Fig. 1 
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Distribution in prong number for (a) stars formed by absorption of the 7-meson 
in a heavy nucleus, and (6) stars formed by absorption of the meson in a 
light nucleus. For light nuclei the results are compared with those of 
Ammiraju and Lederman (1956) for the aborption of mesons in carbon 
and in nitrogen. 


No direct determination was made of the number of events without 
secondary prongs. However, using a mean value of the proportion of 
stars without secondary prongs obtained by Menon et al. (1950), Adelman 
(1952), and Demeur et al. (1956) for a normal emulsion, it is possible 
to determine the total number of negative mesons entering the emulsions. 
Subtraction of the number of events with n>1 then gives the number 
of events with »=0 in the dilute emulsion. The percentages of stars with 
n=0 in light and heavy nuclei have been obtained in this way. 
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In fig. 1(b), along with the results of the present experiment for light 
nuclei, are shown the prong Uistributions obtained for the absorption of 
7 -mesons in carbon and nitrogen by Ammiraju and Lederman (1956) 
using a diffusion cloud chamber. The apparent difference in the per- 
centage of four-prong stars observed by these workers and in the present 


Percentage of stars with m prongs 


Number 
Author of n=0 Al M2 
UT) a 
i | Hola cate riae a iae la 
Menon et al. (1950) 500 — 28 — 8 1h i 13 
De Sabbata et al. 500 — 27 — 6-5 26 | 9-5 9-5 
(1953) 
Demeur e¢ al. (1956) 1038 35 — 3 2 7 10-5 
Present Work 788 TO | 21 2 26 | 10-5 8-5 
+8-5 9 +8 | +10] +3 |+3°5 
| = n=4 C= 
L H L H L 
Menon et al. (1950) BOOM Tes PE ay eee 
De Sabbata et al. 500 12 2-5 4-5 1-5 1-0 
(1953) 
Demeur et al. (1956) 1038 11 2-5 6-5 0-2 1-5 
Present work 788 ll 4 8 —2-5 0-3 
+3 |4+3 [42 |+2-0 |+0°3 |- 
L—Light nucleus H—Heavy nucleus 


experiment is attributable to the presence of a considerable proportion 
of oxygen among the light elements contained in a nuclear emulsion. If, 
as is strongly indicated by the presence of a maximum in the prong dis- 
tribution at »—3, a common consequence of the absorption of a 7~-meson 
by a light nucleus is its complete disintegration into «-particles, a proton 
and neutrons, then oxygen will frequently give a four-prong star. 


ACKNOWLEDGMENT 


The authors are grateful to Dr. H. Muirhead, who originally suggested 
this experiment, for valuable discussions throughout the course of the 


work. 
REFERENCES 


ADELMAN, F. L., 1952, Phys. Rev., 85, 249. 

Ammtirasu, P., and Leperman, L. M., 1956, Nuovo Cim., 2, 283. ; 
Demevr, M., Hutevx, A., and VANDERHAEGHE, G., 1956, Nuovo Cim., 4, 509. 
Dr Sappata, V., Manarest, E., and Purrt, G., 1953, Nuovo Cim., 10, 1704. 
Menon, M. G. K., Murru#ap, H., and Rocuar, O., 1950, Phil. Mag., 41, 583. 


The Nuclear Photoeffect in Light Nuclei; 


By F. C. BARKER 
Research School of Physical Sciences, Australian National University 


[Received February 2, 1957] 


Tue energy of the giant resonance of the nuclear photoeffect is observed 
to be roughly independent of the mass number A for A530, then to 
decrease slowly as A increases (Ferguson ef al. 1954). Several calculations 
of the energy have been made for heavy nuclei, exhibiting the decrease 
for large A (see Levinger 1954). Here we show that the results for light 
nuclei are consistent with the predictions of the nuclear shell model. 

Observationally the most characteristic energy associated with the 
giant resonance is the peak energy H,,. This however is not very 
significant theoretically, and we have calculated instead the harmonic 
mean energy W,, defined by 


[~ #-1¢(B) dE 
Jo 
where o(#) is the total absorption cross section for y-rays of energy F. 
It is difficult to obtain precise values of W, from the observed cross 
sections, as usually only one or other of the partial cross sections, such 
as o,,, has been measured. The estimates of W, which we use for 
comparison with the calculations are given in the last column of the 
table, and are admittedly very rough. This is particularly the case for 
those based on the cross sections of Ferguson et al. (1954), which are 
uncertain at energies above the giant resonance and most likely lead to 
values of W,, which are too small. 
If it is assumed that only electric dipole absorption is operative, then 
(Levinger and Bethe 1950) 
ates (Y,'[DiHD— Di DH]\¥,) i 
H (Pt DIDY,) . : . : . . ( ) 
Here YW, is the ground state wave function, H is the nuclear Hamiltonian 
and D is the electric Ue operator 


v= Stra} $3) 


i=1 
Levinger and Kent (1954) and Levinger (1955) have used shell model 
wave functions in (1) to calculate W,, for various nuclei. An alternative 
formula for Wy is 
(UDHDY,) (PtH) 
(EDDY. (Ee. 
which is identical with (1) if Y is the exact wave function of the ground 


Wy= 
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state corresponding to the Hamiltonian H. For approximate wave 
functions Y, the two formulae give different results. W,, given by (2) 
can be interpreted as the difference H,—H, of the expectation values of 
the Hamiltonian in the state Y,= DY, and in the ground state. 
Only eigenstates of H which contain appreciable amounts of ¥, will 
_be excited by El absorption, and these states will have energies near E 1 
(see Wilkinson 1955). 

For Y, we have taken shell model wave functions belonging to the 
lowest configuration and in L—S coupling, built wp from wave functions 
for single particles in an oscillator potential well with parameter a. 
These wave functions are given in the table, for the various light nuclei 
we have considered. Although they may not represent the actual ground 
states very well, they should be satisfactory for our calculation which 
does not appear to be sensitive to the choice of WY). The nuclear 
Hamiltonian has been assumed to contain a general charge-independent 
two-body central interaction between nucleons : 


Fr A ; 4 
Hf==— Mt Vei+ » [a +a,(o;.0;)+4,(7;.7;) 


4<j=1 
+4y,(0; . 95)(7; - T;)|I (144). 
For all of our calculations we have used a Yukawa potential 
J (r)=V exp (—«r)/Kr. 

If the values of W,, given by (1), (2) are called Wy,, Wy. respectively, 
then we note that matrix elements of only one-particle and two-particle 
operators are required for W,,,;, whereas W,,. involves matrix elements 
of one-, two-, three- and four-particle operators. To calculate Wy, we 
therefore use 

Wie=Wiit 4Wy 
where 
> (Po DUT a) AY) 
4 Ww = nV : 
(¥,D'DY,) 
The orthonormal states Y, include all those which are coupled to Y, 
by both of the operators D'D and H, apart from Y, itself. Because 
oscillator wave functions are used, there is only a finite manageable 
number of such states ¥,. Only one- and two-particle matrix elements 
are required in 4W,,. Thus W,,, and 4W,, can be calculated by standard 
fractional parentage techniques (Elliott 1953). Then 
W a= 77/Ma?-aV 4; 
AW ,=—h?/2Ma?+BV >. 
Tt is found that both « and f are closely proportional to («a)~*, so that 
W y= f?|Ma*(1+-,5), 
W y9=h?/Ma*($+y 2b), 
where b= MJ, /h?«?, and y,=(«a)’« and y,=(«a)?(x+f) are more or less 
independent of xa. Also b does not depend sensitively on « if the 
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internucleon potential is assumed to be the same as that acting between 
neutron and proton in the deuteron, since then (Rosenfeld 1948) 
b=1-680-+2-266n—0-24672+0-14473—..., 
where »=(ME,/f?«?)!?, HE, being the deuteron Pindeay energy. For 
reasonable values of x, ann gives b close to 2-3, which value we have 
used. The calculations have been carried out both for a Rosenfeld 
potential : 
aj=0, dg=0, 4,=0-100, — a,,==0-233, 
and for a Serber potential : 
A=—0'30, @,==0:05, &==0-19, 62 —0-10; 
The value of a has been chosen to make the mean square radius of the 
nucleus equal to 3(r)41/*)?, leading to the values of a/ry given in the 
table. These distances have been measured from the centre of mass 
of the nucleus rather than from the well centre. For the moment we 
assume that r, is the same for all the nuclei considered, and choose the 
value of 7, to make W,, for 1°O fit the experimental value, taken as 
25 Mev. These values of 7) are given in the table below the columns 
giving the values of Wy, and Wyp. 

A further calculation has been made for !2C, including a spin—orbit 
interaction in the Hamiltonian and using intermediate coupling for the 
ground state wave function. It is easily seen, however, that this 
spin-orbit interaction gives no direct contribution to either W,,, (Frankel 
1955) or 4W,. The ground state wave function of “C obtained by 
fitting certain properties of °C and ®N (Lane 1953) is 

YW, —0-829V([44]000) —0-497¥([431]011) +.0-107 ¥([422]000) 

— 0-208 ([422]022) — 0-103¥([332]011), 
all states belonging to the configuration Istlp*. We have approximated 
this by 
Y=0-83¥([44]000) —0-56¥([431]011). 

There are no contributions to W,, from interference between these states 
(such contributions come only from states with the same S and L values 
and even then are very small), so the values of W,, obtained with the 
new Y, are not very different from those obtained previously, being of 
the SHEE of 1 Mev less. 

From the table it is seen that, apart from the values of 7) required, 
the Rosenfeld and Serber potentials give similar results, and the same 
is true of the formulae (1) and (2) for the nuclei with A=10. For these 
nuclei the agreement with the experimental results is Seaeenanie good. 
The discrepancy in the case of °Be may be removed by choosing a larger 
value of ry for this nucleus, of order 1:5x10-%cem; this change is 
supported by experimental evidence from electron scattering and mirror 
nuclei (Hofstadter 1956) which indicates that ry is appreciably greater 
for Li and Be than for He, C or O. For 4He, W,,, seems to agree better 
with the experimental value than does Wy5. If, however, Wj, were of 
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the order of 20 mev, it could not be obtained in the way we have used just 
from the (y, p) and (y,n) cross sections for which the thresholds are at 
about 20 Mey. The values of 7) which we obtain are in agreement with 
the values of about 1-3 10-! em obtained for the nuclei from C to Ca 
from electron scattering experiments (Hofstadter 1956). This success 
of the shell model in fitting the experimental harmonic mean energies 
for the photoabsorption cross sections shows that the shell model is 
capable of giving the energy difference between states belonging to 
different configurations. 

A connection between these calculations and the effective mass concept 
of Van Vleck (1935), Bardeen (1937), Bruechner (1955) and others has 
been pointed out by Professor K. J. Le Couteur. The harmonic mean 
energy calculated for nucleons of mass M* moving in an oscillator potential 
well with parameter a is W,=?/M*a?. Thus by comparison with the 
formula for W,,,, the ratio of the effective nucleon mass to the free 
nucleon mass can be taken as 

M*/M=(1+y,6)71. 

This is about the same for all the nuclei considered here, being 0-6 for 
the Rosenfeld potential and 0-7 for the Serber potential. These values 
are comparable with those deduced from the variation of effective well 
depth with nucleon energy in the optical model of medium and heavy 
nuclei (Melkanoff et al. 1956). The corresponding quantity ($+y,b)"+ 
obtained from W,,. decreases as A increases, being 0-8 for +He, 0-6 for 
*Li, 0-5 for 1®O and 0-44 for #°Ca for the Rosenfeld potential, and slightly 
larger for the Serber potential. 
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SUMMARY 


The angular distributions and energy spectra of the protons emitted,. 
through the interaction of neutrons with an energy of 13-2 Mev with 
natural aluminium, iron and rhodium, have been studied using photo- 
graphic emulsions. The experimental results have been compared with 
the predictions of a model which assumes that these protons can be 
emitted by two processes: by direct collisions of the incident neutron 
with a proton in the target nucleus and by the decay of a compound 
nucleus. Evidence for the process of direct interaction has been obtained 
in all cases. 


§ 1. INTRODUCTION 

IN a previous paper from this laboratory (Morrison et al. 1953) on the- 
interaction of high energy protons with the heavy nuclei (silver and 
bromine) of photographic emulsion, evidence was presented which 
strongly suggested that direct interaction was possible inside the nucleus 
between two nucleons of low energy. The observation of the direct 
emission of protons of low energy was confused, however, by the 
considerable number of evaporation protons emitted in the decay of the 
highly excited nuclei produced in that experiment. It was decided that 
a study of the interaction of nucleons of lower energy could provide more 
definite information concerning the probability of direct emission. 

Paul and Clarke (1953) have made a survey of the cross sections for 
reactions of the type ,N*(n, py), ,N* for neutrons with an energy of 
14-5 Mev, and have shown that the cross sections are, in many cases, 
much larger than would be expected on the basis of the statistical theory 
of Weisskopf and Ewing (1940). These results have been interpreted 
by McManus and Sharp (1955) and by Brown and Muirhead (1957) as 
evidence for direct interaction of the incident neutron with the nucleons 
of the nucleus. Observations on the inelastic scattering of protons 
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with energies of 18 Mev (Gugelot 1954) and of 31 Mev (Eisberg and Igo 
1954) have indicated also that direct interaction as opposed to the 
immediate formation of a compound nucleus is of importance at inter- 
mediate energies. 
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Schematic diagram of apparatus. 


This paper describes a further investigation of the interaction of nucleons 
with nuclei at intermediate energies. The distributions in energy and 
angle of the protons emitted following the interaction of 13-2 Mev neutrons 
with aluminium, iron and rhodium have been obtained using photo- 
graphic emulsions. A study of the (n,p) reaction in other elements 
is in progress. 


§ 2. EXPERIMENTAL PROCEDURE 
2.1. Haposure of the Plates 


The source of neutrons was the T(d,n)4He reaction. The neutrons 
were produced by the bombardment of tritium, absorbed in a target of 
zirconium, with deuterons accelerated in the Glasgow H.T. set. Despite 
the intense flux of neutrons available from this reaction, the small cross 
section for the (n, p) reaction severely limits the possible geometry of an 
experiment in which observations are made on both the angular 
distributions and the energy spectrum of the protons. 


from Aluminium, Iron and Rhodium 787 


_ In view of this consideration the scattering chamber shown in fig. 1 
was constructed. The essential features are the following. Only a small 
portion of the tritium target, 7 in. in diameter, was exposed to the beam 
of deuterons :; thus the neutrons could be considered to arise from a point 
source. Collimation of the deuteron beam above the plate holder ensured 
that the deuterons were incident upon the uncovered portion of the tritium 
target. The scattering chamber was insulated from the H.T. set so that 
the integrated current of deuterons incident on the target could be 
measured. 

Six nuclear emulsions (Ilford C2, 4 in. x 1 in. x 400) were mounted 
hexagonally around this target. Foils, 4 in. x 1 in., of the element to be 
investigated were placed in contact with alternate emulsions. These 
foils and the surfaces of the other three emulsions were covered by gold 
strips of size 4in.x1lin.x400u. The gold, which has a very small 
(n, p) cross section was used to shield the surface of the emulsion from 
external sources of protons. 

Successive irradiations of iron, aluminium and rhodium foils were 
carried out in this arrangement. The time of each exposure was controlled 
by the current integrator. 


2.2. Examination of the Photographic Plates 


The photographic plates were examined in a region corresponding to 
neutrons incident at an angle of 30°--2°. For the exposures with foils 
of aluminium and iron, the area of emulsion examined on each plate was 
6mm x0-3mm,; for rhodium the area examined on each plate was 
6mm x0-6mm. The energy of the neutrons incident over this area 
was calculated from the dynamics of the T(d,n)*He reaction to be 
13-2-L0-2 ev. 

' The plates were examined with «45 and x 90 oil immersion objectives 
and x10 eyepieces for tracks of protons entering the surface of the 
emulsion+. Those tracks whose angle of dip was less than 30° in the case 
of the exposure to iron and less than 45° in the case of the aluminium 
and rhodium were selected for measurement. The limitation on dip 
angle ensured that nearly all tracks that started at the surface ended inside 
the emulsion. The position of the scan was such that measurements on 
protons making angles greater than 150° to the direction of the incident 
neutrons could not be made. Scanning was also carried out inside a 
known volume of emulsion for all protons that satisfied the same geo- 
metrical condition as imposed at the surface. 

- The range, the angle in the plane of the emulsion with respect to the 
direction of the incident neutron and the angle of dip were determined 


+ It was assumed that all tracks of unit charge were protons since it was 
impossible to separate protons and deuterons of low energy by grain counting 
in C2 emulsion. The validity of this assumption will be discussed in the 
presentation of the results. 
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for every track selected. The daily variation of the thickness of the 
-processed emulsion was taken into account in the computation of the angle 
of dip and the range. The energies of the protons were determined from 
their ranges in dry C2 emulsion (Rotblat 1951). The angles in space 
of the protons with respect to the direction of the incident neutron were 
calculated from their dips and projected angles. 


§ 3. ANALYSIS 
3.1. Separation of Events 


Not all protons that were observed to start at the surface of the emulsion 
originated from the elements in contact with the emulsion. Protons. 
starting in the emulsion within a distance of 1 or 2 microns from the 
surface are optically unresolvable from those entering the emulsion. 
Protons in this category were divided into two groups : 


(1) recoil protons produced as a result of the collision of the incident 
neutrons of 13-2 Mev energy with hydrogen present in the emulsion, 

(2) protons produced as a result of the interaction of the incident 
neutrons with the nuclei of the emulsion, or as a result of the collision of 
neutrons of lower energy with hydrogen atoms. These will be called 
background protons. 


The recoil protons (1) were easily distinguished by applying the 
relationship 
H,=(13-2-L0-2) cos? @ 


where #,—energy of the recoil proton and 6=spatial angle between 
the direction of the recoil proton and the direction of the incident neutron. 

Allowance was also made for the errors in the measurements of Z,,. 
and @ in the classification of recoil protons. 

The correction for protons orginating from processes outlined in (2) 
was obtained from the measurements made on protons occurring inside 
the emulsion. Since the ratio of background to recoil protons starting 
at the surface and at points inside the emulsion should be the same, 
the number of surface tracks subtracted as ‘ background’ events was. 
made proportional to that for protons identified as recoil tracks. The 
grouping of the subracted tracks in intervals of energy and spatial angle 
was made in the same relative numbers as those observed inside the 
emulsion. 

The remaining surface protons could have arisen only in the materials. 
placed in contact with the photographic emulsion. The effect of protons 
originating in the gold strip was eliminated by considering the energy 
and angle of the protons left (after subtraction of (1) and (2)) at emulsion 
surfaces in contact with the gold strips alone. Before subtraction from 
the data obtained on emulsions exposed to the foils, allowance was made- 
for the form that this spectrum would take after the traversal of the foil. 
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The procedure described above may be conveniently summarized by 
the schematic equation 


protons from foil= protons observed at the surface 
—|(1)+-(2)4 modified protons from gold]. 


The number of tracks entering the various categories is shown in table 1. 


Table 1 
Element : . ; Protons 
(angle of ee Recoils Saree perch from 
dip) j oat a Foil 
Al 
( ey 1600 220 30 70 1280 
e 
(1 re, 1190 200 (935) 50 885 
(10°—-45°) 1140 490 60 90 500 


3.2. The Correction for the Energy Loss of Protons Produced in the Foil 


The thickness of the foils used in these investigations was 6, 13 and 
12 mg/cm? of aluminium, iron and rhodium respectively. However, 
the protons that were accepted for analysis had traversed the foil obliquely 
so that the effective thickness of a foil was greater than the figure quoted. 
A proton of 5 Mev energy that passes through a foil at an angle of dip 
of 30° loses an energy of 0-7, 1-2 and 1-0 Mev in aluminium, iron and 
rhodium respectively. The energy of a proton produced within a foil 
may therefore be considerably modified, particularly at low energy, 
by its passage through the foil to the surface of the emulsion. For this 
reason protons observed with angles of dip of less than 10° were discarded 
because of the large thickness of foil they might have traversed. 

A relaxation method was adopted to correct for the energy loss of 
protons produced in a semi-thick target. This correction was applied 
for protons originating in the aluminium and iron foils. For rhodium 
the observed energy spectrum was simply corrected to half the thickness 
of the foil, since no protons were found to originate with an energy less 
than 4 Mev. 

3.2. The Determination of the Cross Section for the Emission of Protons 

The differential cross section, (do/dw),, for the emission of protons at 
an angle 0, may be related to the number of protons, V ,(@), observed 
at an angle @ by the expression 

N ,(0)=F (do/dw) pd. 2, nV 
where /=total flux of neutrons per cm? at the foil, d2,—solid angle 
available for acceptance of a proton emitted at an angle 6, n=number 
of atoms per cm? of the foil, and V=volume of the portion of the foil 
under examination. 
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The total flux of neutrons incident on the foils in each exposure was 
determined by counting the number of recoil protons observed within 
a known volume of emulsion. The efficiency for detection of recoil 
tracks was estimated to be 80--10% by examination of the same volume 
of emulsion by different observers. The consequent error in the deter- 
mination of the neutron flux, however, was only taken into account in 
the calculation of the total cross section for emission of protons, since the 
error was systematic and would not affect the relative differential cross 
sections.. 

The solid angle dQ2, was restricted by the geometry of the exposure 
and the limitations imposed upon the angle of dip of the tracks selected 
for measurements. This solid angle was calculated and taken into account 
in the determination of the differential cross sections. 


§ 4. RESULTS 

The energy spectra, over the range from 0° to 140°, of the protons 
emitted from aluminium, iron and rhodium are shown in figs. 2, 3 and 4. 
The corresponding angular distributions for different ranges of energy 
are shown in figs. 5, 6 and 7. The total cross section for the emission 
of protons from each element was obtained by integration of the angular 
distributions extrapolated to 180°. The results are given in table 2 and 
compared with those obtained by other workers. 


Table 2 
Cross section in millibarns 
Nuclear emulsions Activation methods 
Klement 
P aA 
Present nee hae os , 
results : corrected Glen " Forbes 
to n, py nana 
27 Al 140 +30 140 55-15 522510 Sr see 
o6Her 120+30 190 120+30 Sy 124-12 
103 Rh 2bar6 ees ¥ 


+ A target of natural iron was employed in this experiment. To obtain 
the value quoted for °*Fe a small background (~30 mb) due to the presence of 


54 y . “ 
4 He (Allan, private communication) has been substracted from the observed 
cross section, 


The nuclear reactions leading to the production of protons that can 
be induced by neutrons with an energy of 13-2 mev can be represented 
schematically as follows : 


zNA-+n > 7 ,N4+p+y (a) 
yN41+ pn | (0) 
2 NALtntp. @ 
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Fig. 3 
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In the present experiment the sum of the contribution from all three 


reactions is observed. The processes (n, pp) and (n, d) are also 
energetically possible, but, due to the effects of the large negative Q-values. 


Fig. 4 
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Energy spectra of protons from rhodium. 


Fig. 5 
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Angular distributions of protons from aluminium. 


and the coulomb barrier are expected to make a negligible contribution. 
to the observed data. If the reasonable assumption is made that the 
emission of a neutron normally accompanies that of a proton when the 
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energy of the latter is less than 13-2+-q (where ¢ represents the Q-value 
for the processes (b) and (c)) an estimate of the cross section for the 
process (a) can be made by integrating the observed spectrum for energies 
greater than this limit. The cross sections derived in this way are given 
in column 4 of table 2 and compare favourably with those obtained by 
measurements of the activity of the residual nucleus ,_,N* (Paul and 
Clarke 1953, Forbes 1952). 
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Angular distributions of protons from iron. 
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Angular distributions of protons from rhodium. 


In the investigations of Allan (private communication), the angle of 
observation for protons emitted in the three reactions (a), (6) and (c) 
was 34°-20°. When account is taken of the forward peak in the angular 
distribution found in the present experiment the cross sections obtained 
by Allan would be reduced by about 20%. 
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§ 5, Discussion or RESULTS 


The angular distributions at all energies show preferential emission 
of protons in the forward hemisphere. A similar trend in the angular 
distribution of the high energy protons emitted when °*Zn is bombarded 
by neutrons with an energy of 14 Mev has been reported recently (Rosen, 
private communication). A preferential emission of this nature would 
be expected if some protons were emitted following the direct collision 
of the incident neutron with a proton inside the nucleus. The direct 
emission of protons should be enhanced relative to the emission from 
the decay of a compound nucleus for protons of high energy and nuclei 
of high atomic number. It can be seen that, experimentally, the lack 
of symmetry about the position of 90° is particularly evident for protons 
of high energy, and that it is greatest for rhodium. 

It is of interest to compare the experimental results with the predictions 
of a simple model which takes into account both the direct emission and 
that from the decay of the compound nucleus (Brown and Muirhead 1957). 
This model assumes that nucleon-nucleon interactions can occur through- 
out the nuclear volume. 

The magnitudes of the cross sections for the emission of protons in the 
reactions, ,N4(n, py),_,N* and ,N*(n, pn), ,N4~! have been calculated 
using this model. In table 3 the calculated cross sections, o,(dir), for 
the emission of protons by direct collision, and o,(comp) for the emission 
of protons by the decay of an excited nucleus are given. 


Table 3 


Calculated cross section 


aD ee Experimental 
Klement a ie os ie cross section 
og ,(dir) c,,(comp) Total mh barte 
27 A] 80 70 150 140 +30 
56He 25 70 95 120+30 
W3Rh 20 5 25 2545 


The calculated energy and angular distributions are shown in figs. 2-7 ; 
the broken line represents the cross section for the emission of protons 
from the decay of the compound nucleus, while the full curve is the sum 
of this cross section and that for the emission of protons by direct collision. 

The apparent difference in the calculated and experimental energy 
spectra at low energies in aluminium can be mainly accounted for: by 
considering the contribution of protons from the (n, np) reaction. The 
difference between the Q-values for the processes (n, np) and (n, 2n) may’ 
result in the favoured emission of protons over a limited region of. the: 
energy spectrum. This is indicated in table 4. °°. iA ord 
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In the case of aluminium, therefore, the contribution is expected to be 
particularly important. In the case of iron and rhodium the effect is 
reduced considerably since there is no favoured emission of protons over 
the region of the observed energy spectrum. 


§ 6. CONCLUSIONS 


~ The present results for the total cross sections for the (n, p) reaction 
are in satisfactory agreement with those of other workers. Moreover, 
evidence for the process of direct interaction has been obtained in all 
cases from the angular distributions of the protons emitted in the (n, p) 
reaction. 

It can also be seen that the application of the model, proposed by two 
of us (Brown and Muirhead Joc. cit.), to the reactions studied leads to 
rough quantitative agreement with the experimental results. However, 
the theoretical model makes the extreme assumption that a nucleus may 
be treated as a Fermi gas, and that the detailed properties of an individual 
nucleus may be ignored. Thus it is not surprising that detailed 
comparisons do not show more than semi-quantitative agreement. 
Further defects of this model are the neglect of conservation laws other 
than those for energy and linear momentum, and the over simplification 
of the problem of transmission at the nuclear surface. 


Table 4 
Q-values in Mev Maximum proton pee ee 
Element energy from Doe ate 
a ea tn, np) in rev proton emission 
n, np n, 2n ae) 2 in Mey 
| 

27 AT —8-3 —13-0 4-9 0-4-7 
=8he —10-2 —11]-2 3-0 0-1-0 
103Rh —6°5 —9-4 6-7 0-2-9 


The possibility of the enhanced production of protons at the nuclear 
surface also should be considered (Austern et al. 1953). Conservation 
of angular momentum, which is taken into account in this theory, could 
then impose more stringent conditions on the angular distributions of 
the outgoing protons than is assumed in the present model of volume 
production. Unfortunately this theory does not yield absolute cross 
sections. 

It can be seen, therefore, that there is a need to improve the theory 
of nuclear reactions in the intermediate energy region. It must also be 
pointed out that considerable refinements are required for the experiments 
performed in this energy range. None of the experiments so far reported 
for examining the (n, p), (n,n) or (p,m) processes could have detected 
satisfactorily the existence of discrete states superimposed on a continuum 
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as suggested by the experiments of Gugelot (1954) on the inelastic 
scattering of protons. Thus considerable improvement in the energy 
resolution should be sought in future experiments in this field. 
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Observations on Extensive Air Showers 
I. Apparatus} 


By T. E. CransHaw and W. GALBRAITH 
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§ 1. InTRoDUCTION 
AN arrangement of Geiger—Muller (G.M.) counters, covering an area of 
approximately 0-5 km?, together with other apparatus, has been set up at 
Culham (alt. 60m) near Harwell, to study various aspects of extensive 
air showers. 

The results of the investigations will be given in a series of papers, and 
this is intended to be an introduction to the series. It contains a general 
description of the apparatus used, and also includes in Appendices details 
of some of the electronic circuits which may be of interest to other workers 
in this field. Specialized apparatus for particular experiments will be 
described in the appropriate papers. 


§ 2. Larric—E ARRANGEMENT OF G.M. CoUNTERS 


As in an earlier Harwell experiment (Cranshaw and Galbraith 1954) 
G.M. counters on a lattice are used to detect the arrival of extensive air 
showers. The arrangement is shown in fig. 1, and fig. 2 gives a block 


Fig. 1 
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diagram of the apparatus. At each of the ninety-one points in fig. 1 
there is a G.M. unit containing three counters, two of area 200 cm? 
(Cintel G.M. 1A4) and a smaller 15 cm? counter (G.M. 444). The distance 
between units is 90 metres, and the area covered is an equilateral triangle 
of area 0:5 km?. The shape of the lattice was dictated by the area 
available, but it is also convenient for calculation. 

Each G.M. unit contains the following electronic circuits. 

(2) A quench circuit for each counter which lowers the E.H.T. on the 
counter below the threshold for a period of 2 msecs after a dischargey. 
Details of this circuit are given in Appendix I. 

(6) A conventional Rossi coincidence circuit for the two large counters 
with a resolving time of 3 psecs. 

(c) Two blocking oscillator circuits for sending pulses ‘a’ and ‘ b ’ from 
the units to the recording centre at X, fig. 1. 

The ‘a’ pulses are triggered by the coincidences between the two 
large counters. The counting rate of a single counter is about 500/min, 
and the rate of coincidences about 3/min. The ‘b’ pulses are triggered by 
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the small counter, and the rate is about 10/min. The blocking oscillator 
delivers a balanced pulse of 40 »secs duration through a transformer into 
unscreened P.V.C.-insulated flex. The pulse amplitude is about 30 volts 
and the impedance of the line is about 120 Q. At X, another transformer 
converts the balanced pulse into a pulse of single polarity with an ampli- 
tude of approximately 50 volts. Pick up between one G.M. unit and 
another, or between blocking oscillators within the same unit due to 
sending the pulses along the long lengths of unscreened flex is negligible. 
eee 


} At first, a shorter quench time of 80 4 secs was used. The occurrence of 
multiple pulses in some counters after a few months’ operation led us to increase 
the quench time. No trouble has been encountered since then. 
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To test this point, a pulse generator was set up near the centre of the 
array, and leads taken to each of sixteen units so that any combination 
of the units could be triggered simultaneously. Several test patterns 
were tried, and were always registered correctly. 

The positive and negative H.T. supplies and the E.H.T. supply for the 
counters are separately derived at each unit, and the power for the whole 
lattice is provided by a 500 cycle generator stabilized at 240-+-1 volts. 
No attempt has been made to control the temperature of the units, 
because the counters quenched with ethyl formate operate satisfactorily 
over the temperature range encountered in this country. 


§ 3. RecorpING APPARATUS 


The basic circuit used to record the showers is shown in fig. 3. There 
is one such circuit for each G.M. unit. In the pulse adding circuit, the 
‘a’ pulses from each unit are arranged to cut off the current through one 
of ninety-one triodes, whose anodes are connected together. (‘The series 
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germanium diode has the same bias on it as the triode, and prevents small 
pulses picked up on the cables from appearing in the added-up pulses.) 
This current is arranged to flow through a common 5 k2 resistance, across 
which it develops a pulse of about ten volts. Discriminators set at suit- 
able bias levels select showers which discharge 2, 3, 5 and 8 or more G.M. 
units within the resolving time of about 40ysecs. Further details of 
these circuits are given in Appendix II. 

The rate at which showers arrive is recorded by photographing every 
hour in solar and sidereal time a bank of Post Office registers driven from 
the discriminators. Two rows of lights photographed at the same time 
give the barometric pressure and the air temperature on the roof of the 
building, from which the dependence of the shower rates on these two 
variables can be determined. 

For some experiments, we need to know not only how many, but also 
which G.M. units have been discharged. A standard automatic punch 
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(Hollerith type) was adapted to display this information on a punched 
card. When a coincidence is observed, the information is stored on 
miniature thyratrons, type CV797 (see fig. 3). These thyratrons have 
two grids, both of which must be near earth potential for a discharge to 
occur. A master pulse derived from the 3-fold discriminatory is applied 
to the first grid of all the thyratrons, and an ‘a’ or ‘ b’ pulse is applied 
to the other. Then after an event in which three or more ‘a’ pulses are 
produced, there will be at least three ‘a’ thyratrons and possibly some 
‘pb’ thyratrons in the conducting state. A relay circuit designed for us 
by R.C.M. Barnes is then brought into operation by the moving contact 
of the Hollerith machine, and selects all the thyratrons in turn. If the 
thyratron is conducting, an appropriate hole is punched in the card. 
When all the thyratrons have been selected, another relay momentarily 
breaks the H.T. lead to the thyratron anodes so that they are ready for 
the next event. Each shower is recorded on a separate card. 


Fig. 4 
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EVENT NUMBER 
PRESSURE 
“TEMPERATURE 


DAY 


Additional information is punched giving the time to the nearest half- 
hour, the number of the events in the day, the barometric pressure, air 
temperature, and sometimes information from other apparatus. This 
hodoscope unit has operated reasonably well, and the punched cards are 
very convenient for examination and analysis. Figure 4 shows two cards 
as examples. The first card shows a typical three-fold shower, in which 
three ‘a’ pulses and one ‘b’ pulse were recorded. The second card 
shows the largest shower so far recorded. 

a 

_t About one-third of two-fold events are due to random coincidences, and it 
is not worth while to record these on the cards. 


Observations on Extensive Air Showers: I 801 


§ 4. CHECK ON THE CORRECT OPERATION OF THE APPARATUS 

One of the main objectives of the experiment requires reliable and 
continuous operation of the apparatus over long periods. Continuous 
checking was therefore carried out in the following ways. 

(a) The total rate of ‘a’ pulses from all ninety-one units is recorded on 
a chart. Under normal conditions the total ‘a’ rate is 250 counts/min. 

(b) The ‘a’ pulses, and every tenth ‘b’ pulse from each of four G.M. 
units, are recorded on multi-pen charts for periods of half an hour. In 
successive half-hours, four different units are selected, and in this way 
the whole lattice array is checked twice a day. The pulses are taken 
from the point T (figs. 2 and 3) and are fed through a test unit which 
rejects all pulses which have less than a certain amplitude and width. 
The selected values are 30 volts and 30 secs. In this way, it is easy to 
identify units in which the pulse shape has deteriorated. 

(c) The pulse adding circuits and the thyratron storage circuits are 
also checked in sequence at every half-hour interval. A test pulse is 
applied (again at the point T) to four ‘a’ channels and the corresponding 
‘“b’s, and also to another four ‘a’ channels alone. A relay circuit is 
arranged to check that the 8-fold discriminator has been triggered, and 
then a test card with an easily recognized pattern is produced indicating 
that twelve thyratron circuits have operated. Thus test (b) checks the 
apparatus up to the point T, and (c) checks beyond T. 


AP PHAN DX «ol. 
Quench Circuit 

It is well known that the performance of Geiger counters can be greatly 
improved by external quenching. In order that the greatest benefit 
should be obtained, it is essential that the voltage should be lowered as 
early as possible in the discharge, which requires that the trigger circuit 
should be sensitive (Elliot 1949). With most designs of trigger circuit, 
this conflicts with interchangeability of valves and components, because 
the valves are being operated on a part of their characteristics which is 
particularly sensitive to manufacturing tolerances, and to ageing (Cooke- 
Yarborough e¢ al. 1949). In the present experiment, these difficulties 
have been largely overcome by making use of the properties of the diode 
in the negative field region, as discussed by Kandiah (1954). 

The circuit is shown in fig. 5. It will be seen that both halves of the 
double triode are on the working part of their characteristics, some 
asymmetry being introduced in the cathode circuit only to increase the 
size of quench pulse. The loop gain in the circuit is kept below unity by 
the impedance of the diode. The value of this impedance, and there- 
fore the pulse required to trigger the circuit, can be controlled by the 
_current through the diode, i.e. by the choice of Rk. With F fixed, C can 
be chosen to give the required quench time. In our case R=1M, 
C=0-003 nF. The circuit is extremely sensitive, but to a great extent 
independent of valve characteristics. 
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Fig. 5 


280V 


A-PP EON DEX] It 
Pulse Adding Circuits 


The pulse adding circuit is shown in fig. 6. The total current through 
the triodes is 180ma. This current is carried by the 500 2 resistor. 
The anodes are also connected via a condenser to the cathode of a CV2127, 
which is passing 30 ma, and is connected in a voltage stabilizer type of 
circuit, i.e. the impedance at its cathode is approximately 1/9,,G where 
Ym is the mutual conductance of the CV2127, with 5 kQ in the anode, and 
G is the gain of the CV138 stage. Thus the impedance is about 10 2, and 
any current changes in the cathode circuit change the cathode potential 
by very small amounts. Therefore, when the current through the 
triodes is changed by cutting off some of them by ‘a’ pulses, about 2% 


Fig. 6 
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of the current flows through the 500 Q anode load, the remainder flows in 
the cathode circuit of the CV2127, and then develops a pulse across the 
5 kQ resistor in its anode. 

The current in the CV455’s is fixed to within about 1°% by the negative 
line voltage H,,, and the wire wound 0-1% tolerance resistor Rc. The 
pulse developed at the anode of the CV2127 by cutting off n}CV455’s is 
thus nf#,,,Ra/Re. By deriving the bias for the discriminators from a 
potential divider across #,,,, to earth, the setting of the discriminators is 
essentially independent of the value of H,,,, and only dependent on the 
ratio of wire wound resistors. Consistent discrimination between 7-fold 
and 8-fold coincidences has been achieved for two years’ operation. 
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ABSTRACT 


The rates of arrival of extensive air showers near sea level have been 
recorded over a period of almost two years. The results show no evidence 
for sidereal variations in the rates, and the bearing of this in theories of 
the origin of cosmic radiation is discussed. Evidence is presented of 
statistically significant variations in solar time for showers containing a 
mean number of particles of 310’, and the significance of this is briefly 
discussed. 


§ 1. INTRODUCTION 


EXTENSIVE air showers afford the only means of recording the arrival 
of primary cosmic ray particles of the highest energies. Previous 
experiments in this field have shown that the rays exhibit a remarkable 
degree of isotropy up to energies of at least 10!*ev. Evidence for a 
small but none-the-less persistent variation in sidereal time comes from 
the work of Auger et al. (1955) at lower energies (1014-10! ev), and is 
supported to some extent by the early observations of Citron and Kehler 
(1955). Later observations by Citron (private communication) have not 
supported these results and the sum total of events over a period of two 
years do not show any significant sidereal effect. The amplitude of 
Auger et al.’s effect is only a fraction of 1%, but it exhibits maxima when 
the plane of the galaxy passes the zenith at the latitude of observation. 
If one attributes their variation to the distribution of cosmic ray sources 
about the galactic plane, a more striking effect might be expected to show 
up at southern latitudes where the centre of the galaxy passes more 
directly overhead. Farley and Storey (1954) reported a variation with 
an amplitude of about 1% in the rate of showers of about the same energy 
recorded at Auckland, New Zealand, and the period of the maximum 
agreed satisfactorily with the time of transit of the galactic centre. How- 
ever, further observations in Auckland (Farley, private communication) 
appear to show no significant effect to support the earlier published 
result. 

At higher energies, we reported (Cranshaw and Galbraith 1954) an 
experiment at Harwell which seemed to indicate a variation of about 
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5% at an energy then estimated to be about 5x 1018 ey, although no 
variation greater than 0-6°% could be found at 10! ey, and at higher 
energies there were insufficient counts to reveal any variation less than 
13%. Recent work suggests that the data used to estimate the energy 
in these experiments resulted in an overestimate by a factor of 2 or 3. 
It is clear however that at energies up to 101 ev, any variations that may 
exist are less than 1%, and are too small for present experiments to 
demonstrate them beyond question. The great importance of these 
results to astrophysics has led to experiments at higher energies. 


Table 1 


Mean 3 Semi-diurnal term 

Shower | number Rate | ae a ee 
size of per day 

particles g 


Proby. Proby. 


3-fold | 5-7x10® | 70-00-37 : . 0-84 . 0-53 

Sidereal 5-fold | 1-4x107 | 14:5+0-17 S : 0-51 2° 0:29 
time 8-fold | 2:-7x107 | 3-72+0-09 ‘ 2 0:99 : 0-45 
20-fold 10§ | 0-38-+0-03 c 0:76 ; 0-41 


3-fold | 6-7 x 108 2- 0-005 |} 1- 0-40 
2 7 , ; : : 
Sfold | 272x107 | asabove | 9. : 0-01 || 13. 0-003 

20-fold 188 21. 0-18 0-86 


Evidence for a variation in rate of showers of energy 101!” ev has been 
presented by Crawshaw and Elliot (1956). However, in this experiment, 
the rate showed a marked diurnal variation in solar time, and no significant 
variation in sidereal time. Variations in solar time with a large semi- 
diurnal component have been reported by McCusker and Wilson (1956) in 
the rate of dense air showers (>700 particles/m?) and also of penetrating 
showers. 

This paper gives the results of the observations on air showers made at 
Culham in the period December 1954—October 1956 using the apparatus 
described in paper I. We have found no evidence for any variation in 
sidereal time, but in solar time the results suggest a variation which is 
mainly semi-diurnal, and may have some relation to McCusker’s. 
Variations of rate of high energy events in solar time are so unexpected 
that we feel inclined to treat them with suspicion. Nevertheless, they 
have now been reported in several experiments, and we feel it desirable, 
to report our results in some detail. We also discuss the significance of 
the absence of sidereal variations. 


§ 2. ANALYSIS OF RESULTS 


The number of showers discharging n or more G.M. units, with n=3, 5 
and 8, in hourly intervals throughout the solar and sidereal day, was 
obtained from the photographic records of mechanical counters actuated 
by discriminators biased for the appropriate shower sizes. Following 
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the nomenclature of paper I, an event of type n fold or greater, corre- 
sponds to showers which have discharged the two ‘a’ counters in each 
G.M. unit, n or more units being involved. According to paper III, we 
estimate the mean number of particles in the showers detected to be those 
given in table 1. The results of adding the counts over the whole period 
are shown in. fig. 1. 

Fig. 1 
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The data have been analysed in the following way. We first suppose 
that the counts simply obey a Poisson distribution, i.e. that the probability 
of obtaining a count is independent of time. Two tests of this random 
hypothesis are made. Firstly y? is calculated for the hourly totals about 
the mean value, and the results of this analysis are given in table 1. As 
seen there, the only set of observations which is not in reasonable agree- 
ment with the hypothesis is that for the > 8-fold events in solar time. 
Secondly, we calculate the amplitudes and phases of the harmonic 
components present in the data by calculating the best sine waves through 
24 and 12 points for the diurnal and semi-diurnal data respectively, and 
estimate the probability P that such amplitudes should arise on the 
random hypothesis, using the method of Chapman and Bartels (1951). 
The statistical significance of the result is then given by 


P= exp (—e«?)/4) 


where « is the fractional amplitude of the sine wave and & is the total 
number of counts in the data. The results of this analysis are given in 
table 1. 

An analysis of the hodoscope records enabled us to select large showers 
(n>20). In those cases where the shower fell near the edge of the array, 
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we estimated how many units might have been discharged if the lattice 
had been infinite in extent, and so added such events to the category of 
20-fold or greater. The results of a harmonic analysis of this data is also 
given in table 1. 

It is clear from the results of our analysis that there is no evidence for 
any variation in the shower rates in sidereal time. The statistical 
accuracy of the results is such that, setting P—0-1 in the above equation 
as our criterion of significance, we would not consider the variation as 
significant unless its amplitude exceeded 7% for the > 8-fold events and 
25% for the > 20-fold events. In solar time, however, there appears to 
be a semi-diurnal variation in the rate of showers discharging 8 or more 
units, the amplitude being about 13° and the time of the maxima being 
at 02.30 and 14.30hr. The diurnal component also present has an 
amplitude of 9° with a maximum at 15.30 hrs. 


§ 3. Maanetic Riaipiry oF PRIMARY PARTICLES 


It was poimted out earlier (Cranshaw and Galbraith 1954) that any 
anisotropy present in the primary cosmic rays would have to be con- 
sidered in terms of galactic structure, and the interstellar magnetic 
field. The magnetic rigidity Hp of the primary particles therefore plays 
an important part in their motion, and for a particle of charge Z and 
total energy # moving across a magnetic field of strength H is given by 


Hp=E/300Z gauss cm. 


For our largest showers (>20-fold events) we have estimated (paper LIT) 
that the mean shower size is about 9 x 10’ particles, which corresponds to 
a primary energy H,~5x10! ev, if the initiating particle is a single 
proton. The value of Hp for this particle is then 1-5 x10" gauss cm. 
If, however, we suppose the incoming particle is a nucleus of iron, making 
its first collision at the same height as the proton, the number of shower 
particles per nucleon of the primary is 1-6 10°, and allowing for the 
increase in energy per shower particle as primary energy decreases 
(Greisen 1956), the primary energy per nucleon is then ~1-6 x 10" ev 
making H,~9x 1017 ev. In fact, a heavy nucleus makes its first inter- 
action higher on average than a proton, so that the energy per shower 
particle is higher in showers produced by heavy nuclei than in showers 
produced by protons. Thus even for heavy nuclei, a minimum value of 
Hp is about 1-5 x 1014 gauss cm. 


§ 4. DiscussION 
4.1. Absence of Sidereal Variation 


We wish to discuss this result in the light of theories of the origin of 
cosmic rays. 

Following Morrison et al. (1954) we define 5, the anisotropy, as 
2(Jmax—ZJ min)/(Jmax+J min) where Jmax and Jmin are the maximum and 
minimum values of the flux We can set an upper limit to 5 by choosing 
a reasonable value o P in the probability expression (§2). We have 
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taken P=0-1. Then with the numbers of counts, J, in the experiment, we 
arrive at the values of « and hence 5 given in table 2. These values then 
represent fairly generous upper limits to the anisotropy which can be 
supposed to exist undiscovered by the experiment. 

Following Cocconi (1956) we also show in this table the total path 
length, Ltota, which protons of these energies would have if they diffused 
for a sufficient time in a spherical galaxy to acquire the isotropy determined 
by 8, and the estimated maximum path length Lmax before escape from 
the galaxy. As can be seen, the path length requirements to account for 
the observed isotropy are about equal to or exceed the diffusion path 
lengths in the galaxy. This means that even in a spherical galaxy 
containing diffusing elements, protons are hardly retained within the 
galaxy for a time long enough to render them sufficiently isotropic at the 
earth. This argument is even more forceful when applied to a flatter 
disc-like galaxy. 


Table 2 


Ey 


ev € =0-56 Lota aes 
for primary protons ai cm 
ox 10% 0-015 1-5x 1024 | 2x 1024 
10” 0:07 3 1023 | 6x 1023 


5x10" 0:25 10 10% 


A similar conclusion may be reached by considering the motion of 
particles within a spiral arm of the galaxy in a magnetic field of 
5x 10-® gauss. We have shown in § 3 that the radius of curvature of a 
proton of 5 10-1” ev (the energy of the largest showers, >20-fold events) 
would be 3x 10°° cm, which is approximately equal to the radius of the 
spiral arm containing the solar system. Thus it is difficult to see how 
primary protons of this energy could be accelerated and then be retained 
for long enough to acquire isotropy. 

The results of the present experiment therefore favour the suggestion 
put forward recently by Cocconi (1956) and Farley (1956) that the highest 
energy particles, if protons, have an extra-galactic origin. The situation 
for galactic origin is only slightly easier if the incoming particles are iron 
nuclei, for as we have shown, their radii of curvature would be smaller 
than that of protons by less than an order of magnitude. 


4.2. Huidence for Solar Effects 


We wish now to discuss the evidence for solar variations in showers 
which discharge 8 or more G.M. units and make some comparisons with 
other recent experimental work. 

The data contained in fig. 1 have been divided into periods of about 
1 month and the harmonic components calculated. These are shown in 
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fig. 2 on a harmonic dial, each dot representing the amplitude and phase 
of the wave for each month. The radii of the circles are chosen so that 
if the scatter of the points is due to statistical fluctuations only, 39% of 
the points would be expected to lie within the inner circle, and 90% within 
the outer circle. 

It will be seen that the amplitude and phase of the diurnal component 
(fig. 2(a)) varies considerably from month to month, four of the seventeen 
points lying outside the 90% circle. A similar erratic variation was 
observed by Crawshaw and Elliot. 


Fig. 2 
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Approximate dates, 1955: a March, b April, c June, d August, e September, 
f October, g December ; 1956: h January, i February, j March, k April, 
1 May, n July, o August, p September, q October. 


On the other hand, in our experiment the semi-diurnal vectors (fig. 2(b)) 
show a variation which may be reasonably regarded as a statistical 
fluctuation from the resultant vector M, the amplitude of which is 13%. 
This is larger than in Crawshaw and Elliot’s experiment, where the semi- 
diurnal term is about 3%. 

McCusker and Wilson have considered the possibility that the 
variations observed in their experiments may be connected with the semi- 
diurnal pressure oscillations in the atmosphere. The times of the maxima 
in the rate of high density showers agree well with the pressure maxima, 
and the rate of local penetrating showers appears to be in antiphase with 
these maxima. In the present experiment, the times of mimimum rate 
coincide with the times of pressure maxima. 

McCusker and Wilson point out that the variation in rate cannot be 
due simply to variations in the amount of air above the apparatus, since 
they obtain the usual value for the barometer coefficient when short 
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term variations are considered. The same argument applies here. To 
obtain a significant change of rate at sea level from a change in the 
distribution of the air above the apparatus appears to require primaries 
with a very large interaction cross section, larger even than that of heavy 
nuclei. Alfven (1954) has suggested that the highest energy primaries 
might be relativistic dust grains. The behaviour of such particles on 
entering the atmosphere cannot be predicted with certainty, but the 
lateral structure of air showers does not favour this hypothesis (Herlofson 


1956, and paper LV). 
We have as yet no explanation to offer for the variations in solar time. 
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ABSTRACT 


Specific heat measurements between 1-5°K and 20°K, and magnetic 
susceptibility measurements between 1-5°K and 290°K have been made on 
four alloys containing respectively 78-99%, 37-6%, 27:0% and 5:03% by 
weight of cerium in lanthanum. Comparison is made with similar 
previous measurements on pure cerium. The single anomaly in the 
specific heat of pure cerium at 12-3°K splits into two in the alloys, and as 
the dilution is increased these anomalies move to lower temperatures, 
The lower of the two only is accompanied by a change in the magnetic 
properties, a maximum occurring in the susceptibility of the most con- 
centrated alloy and temperature-independent paramagnetism setting in 
suddenly in the next two. In the most dilute alloy the anomalies in both 
the specific heat and in the susceptibility appear to have moved to 
temperatures below 1:5°K. Above 80°K the magnetic susceptibilities of all 
the alloys obeyed a Curie-Weiss law, y=C/(7’—0). The Curie constant 
C varied smoothly from 56-1 x 10-* e.m.u./g Ce for the most concentrated 
alloy (cf. 55-3 x 10-4e.m.u./g for pure cerium) to about 140 x 10-*e.m.u./g Ce 
for the most dilute alloy. 

As in pure cerium, anti-ferromagnetism is believed to be the explanation 
of the anomalous magnetic behaviour of the alloys. However, the 
double specific heat anomalies in the latter suggest that the ordering of 
the 4f electronic orbitals takes place in two stages, while the temperature- 
independent paramagnetism of the more dilute alloys may be the con- 
sequence of freedom of the direction of magnetic ordering to orient itself 
perpendicular to the applied field. The increase in the Curie constant 
suggests that the population of the 4f shell increases, with dilution of the 
alloys, to 2 or 3 electrons per atom of cerium. 


§ 1. INTRODUCTION 
PREVIOUS measurements on cerium (Parkinson et al. 1951, Parkinson and 
Roberts 1957, Lock 1957, referred to hereafter as PSS, PR and L 
respectively), have suggested that its anomalous behaviour below 20°K is 
due to anti-ferromagnetic interactions associated with the 4f shell of 
electrons. In order to investigate the effect of magnetic dilution on these 
interactions, specific heat and magnetic susceptibility measurements have 
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been carried out on a number of alloys of cerium with lanthanum, varying 
in composition from 5% to 80% of the former, at temperatures down to 
1-5°x. Lanthanum was chosen because it forms solid solutions with 
cerium in all proportions (Vogel and Klose 1954), and is only weakly 
paramagnetic (L). 


§2. EXPERIMENTAL DETAILS 
2.1. Specimens 


Four alloy specimens, containing 78-9% (alloy 1), 37-6% (alloy 2), 
27-0% (alloy 3), and 5-03% (alloy 4) cerium by weight, were supplied by 
Johnson, Matthey and Co., Ltd. They were prepared by melting 
together the pure metals, and a portion of each specimen was analysed 
afterwards. The iron content of the specimens was stated to be about 
0:05%. X-ray powder photographs showed that their crystal structures 
were predominantly hexagonal close-packed, like that of lanthanum, 
whereas pure cerium is usually face-centred cubic (Spedding e¢ al. 1956). 

For the susceptibility measurements a small specimen of each alloy, 
weighing 0-4 to 1-1 g according to the concentration of cerium, was cut 
from the main ingot, and the remainder, weighing 35 to 45 g was used 
for the specific heat determinations. 


2.2. Calorimetry 


The specific heats of the alloys were measured from 1-5°K to 20°K in a 
liquefier-cryostat combination (Parkinson 1954), using a thin-walled 
copper calorimeter (PR) employing a constantan heater and a carbon 
thermometer. The specific heat of the empty calorimeter was determined 
in a separate experiment and a small correction was made for the varying 
mass of Wood’s metal used for sealing the lid of the calorimeter. 


2.3. Magnetic Measurements 

The susceptibilities of the alloys were measured from room temperature 
down to 1-5°K, using a Sucksmith balance and cryostat (L). In order to 
correct for ferromagnetic impurities, the measurements were made in 
two different fields (4780 gauss and 2530 gauss) at each temperature, and 
the susceptibility was deduced from the slope of the straight line joining 
these two points on a magnetization curve, plotting magnetic moment 
against applied field. An estimate of the ferromagnetic saturation 
moment can also be obtained in this way. 


§3. REsuLTs 
3.1. Specific Heats 


The heat capacities of alloys 1, 2 and 3 are shown in figs. 1 and 2. The 
atomic heats of pure lanthanum (PSS) and pure cerium (PR) are included 
for comparison. For pure cerium the heat capacity showed a marked 
dependence on the previous thermal history of the sample, and the curve 
plotted is that obtained after the sample had been cooled to 20°k many 
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times. The cerium atoms are then believed to have stabilized with one 
4f electron in each. No dependence on thermal history was observed in 
the present measurements. 

The most interesting feature of the results is that the single specific 
heat anomaly in pure cerium appears to split, in the alloys, into two 
related anomalies. These move to lower temperatures as the proportion 
of cerium is reduced, until in alloy 4 (5-03°% Ce) (not plotted), the observed 
heat capacity is about the same as for pure lanthanum except below 5°K. 


Fig. 3 


Excess heat capacities of alloys per g atom of cerium. — pure cerium ; 
=e Lay A Q, O/ = © 2 / , 
— 78:9% Ce; —-—-37-6% Ce; -— — - —27-0% Ce. 


At this temperature the curve levels off indicating that the anomalies 
probably still exist, but centred below 1-5°K, the lowest temperature at 
which measurements were made. A small irregularity at 3-5°K in the 
curve for alloy 1 may be related to one found at 7°K in pure cerium (PR). 

In order to compare the results for the different alloys, it has been 
assumed that the specific heat of each can be divided into three parts : 
the specific heat of the lattice, the specific heat of the conduction electrons 
and an excess specific heat due to the unfilled 4f shell of the cerium atoms. 


Magnetic Properties of Cerium—Lanthanum Alloys 815 


It is reasonable to assume that the first two of these are about the 
same as in lanthanum (PSS), and therefore the atomic heat of lanthanum 
has been subtracted from the heat capacity of one gram-atom of alloy. 
The difference has then been divided by the atom fraction of cerium in 
the alloy to give the excess heat capacity for a quantity of alloy containing 
one gram-atom of cerium. The results are plotted in fig. 3, together 
with the corresponding curve for pure cerium. The temperatures 7’, 
and 7, at which the maxima in this diagram occur are given in the table. 

Since the temperature scale in fig. 3 is logarithmic, the areas under 
the curves are a. direct measure of the entropy associated with the 
anomalies, and these entropy values are of interest in the interpretation 
of the results. For pure cerium the entropy summation was arbitrarily 
terminated at 20°K, a temperature of approximately 1-6 times that of the 
peak of the anomaly. Therefore the summations for alloys 1, 2 and 3 
have been terminated at about 1:67, that is 20°, 10-5° and 6-5°K 
respectively. These temperatures also coincide with the upper minima. 
in fig. 3. The entropies per gram of cerium, S, thus obtained are given 
in the table. 

The rise in the excess heat capacities for alloys 2 and 3 above 10°K may 
imply the existence of a further anomaly at higher temperatures. A 
similar rise was found in pure cerium above 20°K (PSS). 


3.2. Magnetic Susceptibility 

Figures 4 and 5 show values of 1/y for the alloys plotted against 
temperature, where y is the susceptibility per gramme of cerium, 
corrected for the paramagnetism of lanthanum (L). The correction has 
been made on the assumption that y»=fy+(1—f)x, where y;,, is the 
susceptibility per gramme of lanthanum and x, the observed suscepti- 
bility per gramme of an alloy containing a fraction f by weight of cerium. 
For comparison the values of 1/y for pure cerium after 100 coolings (L) 
are also plotted. 

The ferromagnetic saturation moments, M, of the alloys were roughly 
temperature-independent (except for alloy 3, in which a marked drop 
above 200°K was found) and their values are given in the table. Since 
the saturation moment of pure iron is about 200 e.m.u./g (Bates 1948) an 
iron impurity content of 0-05°%, can give a maximum saturation moment 
of 100x10em.u./g. It is unfortunate that in fact the two most 
dilute cerium alloys gave the highest values of this moment, since the 
corrections to x due to this factor then become exceedingly large, 
particularly in the upper part of the temperature range, amounting to 
30% for alloy 4 at room temperature, and 40% for alloy 3 at 170°. For 
alloys 1 and 2 this ferromagnetic correction is less than 3% at all tem- 
peratures. 

Figure 4 shows that the susceptibilities of the alloys above about 80°K 
obey a Curie-Weiss law, y=C/(7—?), where the Curie constant C 
increases steadily with dilution from 55-3 x 10-4 e.m.u./g for pure cerium 
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to roughly 140 x 10-4 e.m.u./g Ce for alloy 4 (see the table). It is true: 
that the values for alloys 3 and 4 are probably rather inaccurate owing to 
uncertainty in the ferromagnetic correction, but it is improbable that this. 
could account for the entire variation of C observed. It is interesting 
that the temperature hysteresis found in the susceptibility of virgin 
specimens of pure cerium (L) is entirely absent in the alloys. 

Figure 5 shows the behaviour of the low temperature magnetic anomaly 
as the cerium concentration is altered. The maximum in the suscepti- 
bility of pure cerium at 12-5°K moves to a lower temperature for alloy 1. 
In alloys 2 and 3 no maximum is found, but a rather sharp change to. 
temperature-independent paramagnetism occurs. Alloy 4 shows no 
anomaly down to 1-5°K, the lowest temperature at which measurements. 
were made. The temperature at which the behaviour of the suscepti- 
bility changes, 7’,, is given in the table, together with values of the: 
effective atomic moments in Bohr magnetons (u,) calculated from the 
Curie constants per gramme of cerium at the higher temperatures. 


Mxio | Cx108  (cal/°x/ | 


Alloy (e.m.u./g) | (e.m.u./g Ce) (Noe g at Ce) 


Pure Ce 2. == 2-5 55-3 
: pK : - 56-1 


1 ¢ Dek ¢ 
wv oe S ee p 59-2 
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$4, Discussion 


A theory of anti-ferromagnetism in anisotropic crystals has been given 
by Néel (1936) and Van Vleck (1941). It shows that there is a general 
tendency for the atomic moments to line up in an antiparallel array 
perpendicular to the direction of the applied field, and when this happens 
an almost ‘temperature-independent susceptibility should occur. This 
tendency is opposed by the anisotropic crystal field which tends to align 
the atomic moments in a preferred direction in the crystal. When the 
magnetic field is applied parallel to such a direction, the susceptibility 
should fall as the temperature is lowered, tending to zero at 0°K. The 
combination of these two effects in polycrystalline materials gives rise 
to the maximum normally associated with anti-ferromagnetism. Experi- 
mental evidence from susceptibility measurements on single crystals is in 
general agreement with this theory (van den Handel e¢ a/, 1952). 

This suggests a possible explanation for the appearance of temperature- 
independent paramagnetism in the more dilute alloys of cerium with 
lanthanum. In pure cerium and in the concentrated alloy 1, the 
crystalline anisotropy is great enough to cause the usual maximum in the 
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susceptibility when anti-ferromagnetism sets in at the temperature 7's. 
In the more dilute alloys however, although anti-ferromagnetic ordering 
still occurs, its direction is free to orient itself perpendicular to the applied 
field, and this gives rise to temperature-independent paramagnetism. 

If the change in behaviour of the susceptibility at the temperature 7’ 
is indeed indicative of the onset of some kind of magnetic ordering we 
should expect an anomaly in the specific heat at this temperature, and 
we see from the table that the temperature 7’, of the lower maximum in 
the specific heat does in fact coincide very closely with T,. The upper 
maximum in the specific heat, however, remains unexplained. There is 
no observable anomaly in the susceptibility coinciding with it, and yet it 
seems probable that it is closely related to the lower maximum, since 7’, 
and 7’, both decrease steadily with dilution. In addition both specific 
heat maxima must be taken into account in order to obtain an entropy 
contribution per atom of cerium comparable with that for pure cerium. 
It therefore appears that, whatever form the ordering of the 4f electrons 
takes, it is in some way divided into two stages in the alloys, whereas in 
pure cerium the stages coincide. Possibly the difference is associated 
with the fact that the crystal structure of the alloys (hexagonal close- 
packed) has a lower symmetry than that of pure cerium (face-centred 
cubic). 

The increase in the effective moment of the cerium atoms in the more 
dilute alloys suggests that there is a tendency for transitions from the 
outer valency bands of electrons (5d and 6s) into the 4f shell to take 
place, thus increasing the population of the latter. The maximum 
moment found is 4:0u, per atom, slightly greater than the theoretical 
value for neodymium (3-62.,), whose moment should be the greatest 
among the first half of the rare-earth group (Van Vleck 1932). This high 
value may be partly due to an over-estimate of the Curie constant of 
alloy 1 owing to errors in the ferromagnetic correction and in the estimate 
of the paramagnetism of the conduction electrons. The measurements 
suggest however that in dilute alloys of cerium with lanthanum the 
population of the 4f shell is increased from 1 to 2 or 3 electrons per atom. 

Such a transfer of electrons would not necessarily increase the entropy 
in the low temperature anomaly. In pure cerium the degeneracy of 
the ground state of the atom (?F;)) should give rise to a total excess 
entropy of # In (2J-+1)=—3-58 cal deg! (g atom)~1, whereas only a part, 
1-27 cal deg~' (g atom)~ (cf. Rln 2=1-38) appears to be associated with 
the anti-ferromagnetic ordering at 12:5°x. The remainder is found 
between 20°K and room temperature (Skochdopole ef al. 1955). A similar 
entropy distribution may therefore be expected in the alloys, and as 
we have seen, alloys 2 and 3 do in fact show an appreciable rise in 
excess specific heat above 10°K. 

The conclusions reached in this section are of a tentative nature. and 


clearly are not a completely satisfactory explanation of all the features 
observed. 
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ABSTRACT 


A propane bubble chamber of effective volume 1 litre has been con- 
structed for experiments on the elastic scattering of low energy (< 40 Mey) 
a+-mesons on hydrogen and carbon. Stereoscopic pairs of pictures are 
taken on 35 mm film with a stero angle of 90°. The design and operation 
of the chamber and associated equipment is described, and examples of 
pictures taken are shown. 


§ 1, INTRODUCTION 


THE mechanism governing the operation of bubble chambers is not fully 
understood but a theory proposed by Glaser (Glaser 1953, 1954) enables 
the approximate operating conditions to be determined for a variety of 
liquids. Since the pioneer work of Glaser in 1952 (Glaser 1952) several 
authors (Glaser 1953, Glaser and Rahm 1955, Leitner et al. 1955, 
Parmentier and Schwemin 1955, Bassi et al. 1956, Glaser and Brown 1956, 
Glaser et al. 1956, Nagle et al. 1956 Pless and Plano 1956) have published 
details of bubble chambers of the two main types in current use, namely 
“clean ’ chambers, constructed entirely of glass and ‘ dirty ’ chambers 
fabricated of metals, gaskets and glass. 

The present paper gives details of the construction and operation of 
a ‘ dirty ’ propane bubble chamber which has been designed to investigate 
the interactions of low energy pions (<40 Mev) with hydrogen and carbon 
using the external pion beams of the Liverpool synchrocyclotron. At 
these energies the range of the recoil protons in propane is short. To 
minimize the risk of missing short recoil tracks during the analysis of the 
pictures two plane projections of the sensitive volume of the chamber at 
90° to each other are photographed using dark field illumination. This 
procedure also makes the computation of the true geometry of the inter- 
actions relatively simple and eliminates corrections for the high index of 
refraction of the glass windows and the chamber liquid. 


§ 2. DESCRIPTION OF THE CHAMBER 


The body of the chamber (fig. 1 (a)) is a rectangular block of MG5 
aluminium alloy, 6in.x6in.x12in. with a 3in. diameter cylinder 
removed along the axis of greatest length. Two rectangular windows of 
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}in. armour-plate glass are fastened to adjacent faces of the chamber 
body by flanges of MG5 using Bridgman type seals with ‘ Gaco ’ oil-seal 
rubber gaskets (fig. 1 (b)). Each window is 6in.x3in. and a useful 
volume of approximately 1 litre is photographed. In this region the 
inner diameter of the chamber is increased to 32 in. to prevent the direct 
illumination of the walls. The illumination window is 3 in. in diameter 
and is attached in the same way as the other windows. Opposite this is 
a thin metal window of 2} in. diameter and 0-194 in. thick through which 
the beam of particles enters the chamber. The whole of the inside 
surface of the chamber is highly polished to reduce wall boiling and 
anodised black to provide a good photographic background. When 
fully assembled the chamber was tested with a static pressure of 800 p.s.i. 
which is twice the operating pressure. 

The chamber is pressurized by means of a free-moving aluminium piston 
(Leitner et al. 1955) running in a brass cylinder with a polished chromium- 
plated liner. Sealing is effected by means of two ‘ Gaco ’ piston seals, the 
space between them being connected to the atmosphere. The position of 
the piston is measured by pushing a rod into contact with its outer end, 
and this enables the change in volume during an expansion to be measured. 
The expansion cylinder communicates with the bubble chamber through 
a low impedance ‘ neck’ 44 in. long and 2 in. in diameter. 

The expansion cylinder, safety valve, condenser-type pressure gauge, 
and filling lne are attached to the chamber in the unphotographed 
region to reduce the amount of background boiling appearing on the 
photographs. 

The filling system allows propane gas to pass through a filter and 
condense in the chamber. There are facilities for evacuating the chamber, 
adjusting the amount of propane therein and for measuring the static 
pressure. 

The whole of the chamber and expansion cylinder are immersed in a 
water bath, of capacity ~100 litres, the temperature of which is thermo- 
statically controlled to within -+-0-2°c in the range 56°c to 62°c. 


§ 3. THE EXPANSION SYSTEM 


To operate a bubble chamber of the ‘dirty’ variety the expansion 
system must be capable of performing, in a reproduceable manner, a 
rapid expansion and a rapid recompression a few milliseconds later. 
A rapid expansion is necessary since the minimum pressure obtained during 
an expansion is the net result of a lowering of the pressure, due to 
movement of the piston, and a return to the saturated vapour pressure 
due to boiling on the walls and gaskets. If the bubbles formed are 
permitted to coalesce at the top of the chamber it may be necessary to 
delay recycling for several minutes while these bubbles collapse, but, if 
the recompression is sufficiently rapid most of the bubbles can be made 
to collapse before they reach the top of the chamber (Leitner ef al. 1955). 
Rapid collapse of the bubbles is assisted by pressurizing the chamber 
above the saturated vapour pressure at the operating temperature. 
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Fig. 1 (b) 
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A commercial valve with the required speed and reliability could not be 
found in this country when the work was started. The expansion system 
used (fig. 2) consists essentially of two valves. The first C is a high- 
pressure high-speed three way valve of novel design which is connected to 
the expansion cylinder. This valve is in turn controlled by a low-pressure 
three-way pilot valve. The piston of the high-pressure valve C is held 
sealed on a } in. diameter hole, against the 375 p.s.i. in the expansion 
cylinder, by applying 65p.s.i. to its larger diameter end. When the 
lower pressure is released and falls below that required for sealing, the 
high pressure acts on the full area of the | in. diameter end of the piston 
of valve C. After the first } in. of its stroke this piston imposes a high 
impedance in the high pressure gas inlet. During the second } in. of its 
stroke, when it is moving fastest, it opens two } in. diameter holes to the 
atmosphere releasing the gas acting on the piston of the expansion 


Fig. 2 
VALVE A 
4 


VALVE B 


VALVE C 


EXPANSION 
CYLINDER 


BALLAST 
CYLINDERS 


Expansion system. 


cylinder and thus expanding the chamber. The replacement of the low 
pressure gas in valve C returns the piston of C to its original position and 
thus re-applies the high pressure gas to the expansion cylinder. Sealing 
takes place in this valve before the pressure in the expansion cylinder 
has fully built up. 

The low pressure gas in valve C is controlled by valves A and B. A is 
a commercial solenoid operated valve. On energizing or de-energizing 
this solenoid the gas pressure from one end of valve B is removed and 
pressure is applied to the other end. This causes the piston of valve B to 
move from one end to the other. During its movement this piston 
releases the gas from the low pressure end of valve C and replaces it after 
a short delay. Valve B although constructionally asymmetric is sym- 
metric in operation. By varying the low pressure applied to valves A 
and B the time for valve B to complete its stroke is altered, thus changing 
the delay between expansion and recompression of the chamber. 
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In this way the chamber is made to undergo a complete expansion and 
recompression cycle by either energizing or de-energizing the solenoid of 
valve A. Due to the high impedance of valves A and B about 80 msec 
elapses from the start of operation of valve A to the beginning of the 
expansion stroke in the chamber. Most of this delay is due to the high 
impedance in A but the delay is constant to less than } msec. T'wo gas 
pressures are required, one (~65 p.s.i.) for valves A and B and the other 
(375 p.s.i.) for valve C; both of these are supplied by one compressor 
which feeds two automatically adjusting reducing valves. To avoid 
variations in pressure due to gas flow through long supply lines, ballast 
cylinders are placed close to the expansion valve system, the connections 
to the latter being made by short pipes of wide bore. 


§ 4. PRESSURE MEASUREMENT 


A capacitance transducer gauge with a frequency response of 25 ke/s 
and a linear capacitance change of 20 pr over the range of pressures 0 to 
720 p.s.i. is fitted into the wall of the chamber (fig. 1 (a)). The variations 
in capacitance of this gauge due to changes in the pressure in the chamber 
are used to modulate the frequency of a Franklin oscillator (fig. 3). The 
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resulting frequency changes are transformed into voltage changes by 
means of a ratio detector (Sturley 1955) and after amplification these 
voltage variations are fed on to the Y plates of a specially designed 
oscilloscope. To ensure that the full pressure cycle takes place with the 
ratio detector working in its linear region the frequency of the oscillator 
is set at 10-7 Mc/s when the chamber pressure is 230 p.s.i. 

The oscilloscope used has a triggered time base of 200 msec length and 
the tube has a persistence of 10 seconds. The brightness can be modulated 
by a multivibrator of period 0-75 msec thus providing a time calibration 
on the trace. It also has facilities for mixing sharply defined timing 
pulses with the output from the pressure measuring circuit. : 
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The trace is triggered at the same time as the expansion valve begins 
to operate and after a fixed interval (~80 msec) the variation of the 
chamber pressure with time is displayed. Figure 4 (Plate) shows a 
photograph of the trace obtained. The full cycle of expansion and 
recompression takes place in <30 msec, the lowest pressure of about 
180 p.s.i. is attained in 3 msec and the chamber is sensitive for less than 
2msec. It is seen from the photograph that much of the re-pressurization 
is due to boiling in the chamber and occurs before the high pressure gas is 
re-admitted into the expansion cylinder. 

The output voltage from the pressure measuring circuit is also used 
to measure the static pressure in the chamber by applying it to a volt- 
meter calibrated in p.s.i. This measurement is important since the 
chamber is isolated from the pressure gauge in the filling system during 
operation. 

Great care is taken to ensure that the dynamic pressure changes in the 
chamber are accurately represented by the oscilloscope trace. A know- 
ledge of the minimum pressure attained during an expansion at a given 
temperature is extremely useful since the minimum pressure required for 
sensitivity at a given temperature can be calculated approximately from 
Glaser’s formulae (Glaser 1953). It is found that the minimum pressure 
obtained with the given expansion system, rises more slowly with tempera- 
ture than that given by Glaser’s formula. By comparing the variation 
in the minimum pressure obtained as the temperature is raised with that 
required, it is possible to predict approximately the minimum operating 
temperature for a given instrument. 

An accurate knowledge of the dynamic pressure in the chamber also 
enables the effect of varying the speed of expansion to be assessed 
quantitively, and is of great value during the initial testing of the bubble 
chamber. 


§ 5. THe TiImiInG SEQUENCE 

During normal operation the 156 in. synchrocyclotron provides a beam 
of accelerated particles in short bursts (100 to 400 usec long) repeated at 
10 msec intervals. For use with the bubble chamber the machine can be 
operated to give a single burst of particles. The sensitive time of the 
bubble chamber is not more than 2 msec and to ensure that the beam of 
particles traverses the chamber during this interval the operation of the 
chamber is synchronized with the modulation frequency of the synchro- 
cyclotron radio frequency system. 

An optical system attached to the shaft of the rotating condenser of 
the synchro-cyclotron generates a positive pulse for each modulation cycle 
of the radio-frequency. One of these synchronizing pulses (fig. 5) is 
arbitrarily selected by the bubble chamber control unit and successive 
pulses are counted using a cold cathode decatron tube (CV 2223). After a 
number of pulses x, not greater than 10, have been counted, the (n+-1)th 
pulse is used to initiate a sequence of operations which causes the synchro- 
cyclotron to accelerate a single pulse of particles after a delay of 13 msec. 
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The total time 7' between the selection of the first synchronizing pulse and 
the entry of the particles into the chamber is thus 10(n+-1)+138 msec. 
The time taken for the chamber to attain a sensitive condition after the 
actuation of the expansion system is (7'—t) msec and by a suitable choice 
of n, t can be made not more than 10 msec. The expansion system is 
therefore pulsed ¢ msec after the first synchronizing pulse is selected, t 
being continuously variable from 0 to 10 msec. 

The entry of the beam of particles into the chamber is registered by 
two scintillation counters in coincidence. The coincidence pulse is used 
as a time marker which can be mixed with the output of the electronic 
pressure gauge and appears as a sharp pulse superimposed in the pressure 
time curve. To obtain the correct timing sequence for operation of the 
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chamber the oscilloscope is observed and the two variables n and ¢ adjusted 
so that the time marker pulse appears on the correct region of the pressure 
time curve. The coincidence pulse also triggers the light source after a 
delay which is variable from 0 to 3 msec. 

The complete sequence of operations described above is regularly 
repeated at predetermined intervals variable from 2 to 40 sec. 


§ 6. ILLUMINATION AND PHOTOGRAPHY 


The camera and illumination are arranged to give a 90° stereoscopic 
angle and dark field illumination with the direction of the light at 90° 
to the directions of photography (fig. 1 (b)). A high intensity light source 
is required since only a small fraction of the incident light is scattered 
into the camera lens. The source must also be of small dimensions, since 
the light traversing the chamber must be well collimated to minimize the 
direct illumination of the walls. 

A light pulse of 20 jzsec duration is produced by discharging a condenser 
of capacity 33 uF charged to 4000 volts, through a spark gap in air. This 
spark gap is a type described by Barddécz and Klatsmanyi (1955) and can 
be easily triggered. A condensing lens placed in front of the spark gap 


Propane Bubble Chamber 827 


produces a beam of light which is slightly convergent. To increase the 
intensity of the light a concave mirror is placed behind the spark gap 
and the inner surface of the beam entry window is chromium plated and 
highly polished. 

A camera and mirror system are mounted as shown in fig. 1 (b). The two 
plano-convex lenses are placed so that the camera lens is in their focal 
plane (Campbell and Pless 1956). With this system two plane projections 
at right angles are obtained side by side on the film. Photographs are 
taken with a camera lens of 3 in. focal length working at //8 on Kodak R60, 
35 mm recording film giving a demagnification of 8-5. 

For analysis each pair of pictures is projected on to a translucent screen 
to give images double the real size. Measurements are made on each 
projection separately by viewing against the direction of the light and 
the true measurements are computed. 


§ 7. CHARACTERISTICS OF THE CHAMBER 


In the operation of a bubble chamber there are five principal parameters, 
namely: (1) the temperature, (2) the minimum pressure reached during 
an expansion, (3) the volume change during an expansion, (4) the static 
pressure applied to the propane, (5) the delay between the entry of the 
particles into the chamber and the taking of the photograph. The effect 
of varying these parameters on the operating characteristics has been 
found and will now be discussed. 

Below 56°c the chamber is not sensitive because the minimum pressure 
reached is not far enough below the vapour pressure to produce a sufficient 
degree of superheat. As the temperature is raised to 56°c the chamber 
becomes sensitive to the electron track endings produced when a ®Co 
source is placed near. The volume change needed for this condition to 
be attained is 3°% of the chamber volume, and this is called the expansion 
ratio. As the temperature is raised to 59°c, the chamber records 95 Mev 
a-mesons (B=0-8) with an expansion ratio of 2°% and a minimum 
pressure of 170 p.s.i. At 61°c the chamber is sensitive to minimum 
ionizing particles; the minimum pressure under these conditions is 
180 p.s.i. An increase in the expansion ratio from 2% to 4% does not 
alter the minimum pressure reached but increases the sensitive time. 
As the timing sequence is reproduceable to } msec an expansion ratio of 
2% is sufficient. 

The effect of increasing the static pressure applied to the propane is 
to increase the possible recycling rate by recompressing more rapidly the 
bubbles formed during an expansion. At 61°c the vapour pressure of 
propane is 320 p.s.i. and with a static pressure of 370 p.s.. the chamber 
can be recycled every 5 seconds. 

The time taken for the bubbles to grow to a convenient size for photo- 
graphy is 2 msec, but photographs have been obtained with flash delays of 
300 psecs. With a delay of 2 msec no distortion of the tracks of 95 Mev 
=-mesons is observable, other than that due to multiple coulomb scattering. 
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These results were obtained with a filling of commercial quality pro- 
pane, the impurities being ethane, 2% and sulphur, 0-003%. 

Apart from the parameters discussed the characteristics of the chamber 
depend on the smoothness of the internal surfaces. When initial attempts 
were made to operate the chamber, the sensitive volume was seen to 
contain large ‘plumes’ of bubbles in addition to tracks (fig. 6, Plate). 
This effect is thought to be due to roughness of the chamber walls as a 
result of the machining because the effect is no longer observed since the 
inner surfaces have been highly polished and re-anodized (fig. 7, Plate). 


§ 8. Discussion 


A bubble chamber has- several advantages over a diffusion cloud 
chamber or a Wilson cloud chamber for the proposed experiment on low 
energy pion scattering. The medium has a high density (0-43 g/cm?) and 
a larger number of interactions will be obtained per picture than in a 
cloud chamber. The rate at which pictures can be taken is also greater. 

The bubble chamber does however have several disadvantages. The 
medium contains both hydrogen and carbon (C,H,) and the interactions 
occurring in H and C have to be identified accurately from the kinematics. 
The high stopping power of the medium makes the recoil tracks very 
short and therefore difficult to detect and measure. This latter difficulty 
is partially overcome by using 90° stereo-photography and cumbersome 
corrections for the high refractive index of the glass and propane are 
avoided by using plane projections. 

A hydrogen bubble chamber is an ideal instrument for this experiment 
but is considerably more difficult to construct and operate. 

No magnetic field is applied to the propane chamber since the multiple 
coulomb scattering is large at these low energies and a very large field is 
required (~30 k gauss). 

Recycling rates considerably faster than one every 5 sec are mechani- 
cally possible. The recycling rate used at present is determined by the 
amount of high pressure gas available. 

The instrument has been used for several days in conjunction with the 
156 in. synchrocyclotron. The initial setting up is accomplished by 
observing the trace on the oscilloscope and varying the delay between the 
start of the expansion valve cycle and the entry of the particles into the 
chamber until the coincidence pulse appears on the correct portion of the 
pressure time curve (fig. 4, Plate). Tracks can also be observed in the 
chamber visually. 

The optical resolution achieved with the system used is not high enough 
for accurate bubble counts to be made. This is not however a serious 
disadvantage for the proposed experiment since at low pion energies the 
tracks are essentially solid unless extremely high resolution is used or 
the chamber is operated at low sensitivity (Glaser et al. 1956, Blinov 
et al. 1956). The energy spectrum of the incident z-mesons can be 
obtained by measuring the track lengths in the chamber and the energy 
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of scattering events determined from the kinematics. Bubble counting is 
very tedious and time consuming and would not in any case, be used 
when several thousands of tracks have to be measured. It is likely to 
be of use in the future to assist in the identification of particles from their 
tracks in bubble chambers (Budde et al. 1956). 
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REVIEWS OF BOOKS 


Thermodynamics and Statistical Mechanics. By A. H. Witson. (Cambridge 
University Press, 1957.) Price 50s. 

Tus is an advanced textbook written, as the author says, for theoretical 
physicists. It will be invaluable for them and for the growing number of solid- 
state physicists, physical chemists and others who need a thorough understanding 
of both thermodynamics and statistical mechanics and of their relationship to 
atomic mechanisms in order to understand the problems that they investigate. 
The author develops these subjects and has applied them to selected fields, such 
as electrolytes and electrochemical systems, chemical reactions and superlattices 
in alloys. As a clear and reliable introduction, the book can be thoroughly 
recommended. N. F. M. 


Radiation Handbook. (National Bureau of Standards Handbook 62. US. 
Government Printing Office, Washington, D.C.) Price 40 cents. 


THis small volume published by the International Commission of Radiological 
Units and Measurements contains much information basic to the practice of 
dosimetry as well as the most recent international recommendations. It can 
confidently be commended to the attention of readers of the Philosophical 
Magazine. 7H G. 


[Lhe Editors do not hold themselves responsible for the views 
expressed by their correspondents ] 


R. C. NEWMAN Phila Mag. Ser. 8, Vol.-2,-Pl, 26. 


Vig. 1 Fig. 2 


Silver substrate plus 2 A of lead: silver 6 A of lead on silver: (111) plane, [110] 
[110] azimuth. azimuth. 
Fig. 3 


Silver substrate plus 1-5 4 of thallium : 8 A of tin on silver: silver [110] 
silver [110] azimuth. azimuth, 


R. C. NEWMAN Phil. Mag: Ser. 8) Vol 2, aPle 27; 


Fig. 5 


60 A of tin on silver: silver [211] 40 A of silver on silver: (111) plane, 
azimuth. [211] azimuth. 
Fig. 8 


10 4 of antimony on silver: silver [110] 8O A of antimony on silver : pattern 
azimuth. independent of azimuth setting. 
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Fig. 4 


Al 'B Ct 
A photograph of the pressure time curve displayed on the oscilloscope. 
A. The chamber starts expanding. 
B. Particles enter the chamber. 
C. High pressure is re-applied to the chamber. 


Fig. 6 


photograph of low energy 7+-meson tracks in the bubble chamber, 


A single 
showing the ‘ plumes’ produced by boiling on the walls. 
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A pair of photographs showing 7*-mesons stopping in the bubble chamber. 


